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ABSTRACT 



Various extensions of the standard model of particle physics predict the existence 
of very light bosons, with masses ranging from about 10 -5 eV for the QCD axion down 
r \ , to 10~ 33 eV for ultra-light particles. These particles could be responsible for all or part 

of the cold dark matter (CDM) in the Universe. For such particles to serve as CDM, 
i-C | their phase-space density must be high enough to form a Bose-Einstein condensate 

Ph, (BEC). The fluid-like nature of BEC-CDM dynamics differs from that of standard 

q collisionless CDM, however, so different signature effects on galactic haloes may allow 

$— i \ observations to distinguish them. Standard CDM has problems with galaxy observa- 

■ tions on small scales; cuspy central density profiles of haloes and the overabundance of 

subhaloes seem to conflict with observations of dwarf galaxies. It has been suggested 
that BEC-CDM can overcome these shortcomings for a large range of particle mass m 
and self-interaction coupling strength g. For quantum- coherence to influence structure 
on the scale of galactic haloes of radius R and mass M, either the de-Broglie wavelength 
\Q '. X deB < R, which requires m>m H = l(T 25 (i?/100 kpc)- x / 2 (M/10 12 Mq)" 1 / 2 eV, or 

\f} ' else XdeB *C R but gravity is balanced by self-interaction, which requires m 3> m# and 

g > g H = 2 • 10" 64 (i?/100 kpc)(M/10 12 M©)" 1 eV cm 3 . Here we study the largely- 
neglected effects of angular momentum on BEC haloes. Dimensionless spin parameters 
' ^ — 0-05 are expected from tidal-torquing by large-scale structure formation, just as 

, for standard CDM. Since laboratory BECs develop quantum vortices if rotated rapidly 

enough, we ask whether this amount of angular momentum is sufficient to form vortices 
in BEC haloes, which would affect their structure with potentially observable conse- 
quences. The minimum angular momentum required for a halo to sustain a vortex, 
Lqm, corresponds to H per particle, or HM/m. For A = 0.05, this requires m ^ 9.5m//, 
close enough to the particle mass required to influence structure on galactic scales that 
BEC haloes may be subject to vortex formation. While this is a necessary condition, 
it is not sufficient. To determine if and when quantum vortices will form in BEC halos 
with a given A-value, we study the equilibrium of self-gravitating, rotating, virialized 
BEC haloes which satisfy the Gross-Pitaevskii-Poisson equations, and calculate under 
what conditions vortices are energetically favoured, in two limits: either just enough 
angular momentum for one vortex or a significant excess of angular momentum. For 
A = 0.05, vortex formation is energetically favoured for L/Lqm 5^ 1 as long as both 
m/rriH ^ 9.5 and g/gn ^ 68.0. Hence, vortices are expected for a wide range of BEC 
parameters. However, vortices cannot form for vanishing self-interaction (i.e. when 
XdeB ;$ R), and a range of particle parameters also remains even for BEC haloes 
supported by self-interaction, for which vortices will not form. Such BEC haloes can 
be modelled by compressible, (n — 1) -poly tropic, irrotational Riemann-S ellipsoids. 

Key words: cosmology: theory - dark matter - galaxies: haloes - galaxies: kinematics 
and dynamics - methods: analytical 
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1 INTRODUCTION 

In the last decades, astronomical observations have provided 
a range of supporting evidence that about 23 % of the en- 
ergy density of the Universe is comprised of non-baryonic 
dark matter. Its particle nature is still unknown, but all 
candidates find their justification in extensions of the stan- 
dard model of particle physics. Based on observations, it 
is now believed that dark matter is weakly-interacting and 
non-relativistic, hence behaving like a cold gas, and gener- 
ally termed cold dark matter (CDM). 

The standard form of CDM is often supposed to be a 
relic of the Big Bang in the form of weakly-interacting, mas- 
sive particles ( WIMPs) , in particular, the lightest supersym- 
metric particles, the most popular of which is the neutralino, 
with particle mass of the order of 100 GeV. Efforts are un- 
derway to measure the presence of those particles, but no 
direct detection has yet been reported. On the other hand, 
as far as the dynamics of galaxies and large-scale structure is 
concerned, standard CDM is modelled as a cold, collisionless 
gas with vanishing pressure. While many observational prop- 
erties of galaxies can be reproduced, some crucial issues are 
subject to controversy. Simulations of standard CDM struc- 
ture formation predict a universal halo density profile, which 
has a cusp going like r _1 in the center. This seems to be 
in conflict with the observational properties of low-surface 
brightness (LSB) and dwarf galaxies. Moreover, standard 
CDM predicts the hierarchical merger of smaller structures 
into larger structures over time, with an overly abundant 
population of subhaloes inside a halo like the Local Group 
for which there does not seem to be an observational coun- 
terpart. Given the above shortcomings, and the null results 
of direct detection experiments to date, we are still free to 
consider other candidates for the CDM paradigm. 

Another major candidate for dark matter is the QCD 
axion, the pseudo-Nambu-Goldstone boson of the Peccei- 
Quinn phase transition, proposed as a dynamical solution 
to the CP-problem of the strong interactions. For the ax- 
ion to be CDM, it has to be very light, m ~ 10 -5 eV, 
and direct detection experiments rely on its ele ctrody- 
namic al coupling, employing the h aloscope idea o f ISikiyid 
l|l983h in the ADMX experiment (|Asztalos et alj l^Old ) 1 ). 
and axion-photon oscillations in the presence of high mag- 
netic fields using l asers (so-cal l ed li ght-shining-through-wall 
experiments), e.g. IChou et al.l (2008). 

In addition to the QCD axion, multidimensional 
cosmological and string theories generically predict 
the existence of even much lighter bosonic particles 
down to masses of th e order of 10~ 33 eV (see e.g. 
iHorava fc Witter] (jl996t ): iGiinther fc Zhukl jl997r i) which 
could for m all or part of the dark matteij, e .g. ICarrollI 
Il998t): [Arkani-Hamed. Dimopoulos fc Dvalil l|l999l) 



l Arvanitakietah l For such particles t o se~as 



CDM, their phase-space density must be high enough to 
form a Bose-Einstein condensate (BEC), i.e. a macroscopic 
occupancy of the many-body ground state. Generally, 



1 While those particles all the way up to the QCD axion are 
often summarized as 'axion-like' particles, we will avoid this term 
here, since the definition of what is 'axion-like' may sometimes 
encompass particles in the keV mass range. Instead, we will refer 
to them as ultra-light bosons. 



for a system to undergo Bose-Einstein condensation, the 
phase-space density n\ deB must exceed a number of 
order one, where n is the number density and XdeB is the 
de-Broglie wavelength of the particles. While for axions 
n ^deB 3> If ultra- light bosons even fulfill n\ 3 drR S s > 1. I t 
has been shown in mod els bv ISilverman fc Malletl (2002), 
lUrena-Lopej (|2009l ) and iLundgren at alj (|2010h that ultra- 
light particles are able to undergo a phase transition in the 
early universe to a BEC at a temperature which is well 
above that of the time of recombination. These bosons are 
thereby able to form a non-relativistic, quantum-degenerate 
gas, in contrast to the (fermionic) neutrinos which remain 
relativistic. While numerous searches for the QCD axion are 
currently pursued, the ultra-light bosons seem to be out of 
reach for direct detection. However, low-temperature ultra- 
precision experiments promise new possibilities to search 
directly even for these ul t ra-ligh t particles, for a review 
see e.g. Ijackel fc Ringwaldl (|2010h . An interesting, but very 
challenging experimental proposal for detecting such parti- 
cles, by exploiting their sup posed coupling to gluons, ha s 
been presented recently by iGraham fc Raiendranl (|20 1 lh . 
Nevertheless, this kind of dark matter will best be traced in 
the near future through the signature of its dynamical differ- 
ences from standard CDM, an aspect of which we will study 
in this paper. A full account of all the existing literature is 
beyond the scope of this paper, but earlier investigations of 
mo dels for light bos o nic da rk matter include e.g. the works 




Ipser fc Sikiviel (|l983l); iBaldeschi. Gelmini fc RuffinS 



Khlopov. Malomed fc Zeldovich 



iPress. Rvden fc Spergell 



Membrado. Pacheco fc Sanudel (Il989h 



Widrow fc Kaiserl 
Lee fc Kohl ll 19961): 



ISchunckl jl998l); _ 
ISchunck. Fuchs fc Mielkel l|2006h and lMatos fc Urena-Lopej 
(2001) to name a few. Most of this literature has been 
restricted to free particles without self-interaction. 

We will see that BEC cold dark matter (henceforth 
also BEC-CDM) obeys quantum-mechanical fluid equations. 
Therefore, small-scale structure can be very different for 
BEC dark matter from that of collisionless standard CDM. 
The behaviour of self -interacting BEC-CDM as a super- 
fluid makes possible entirely new phenomena like the for- 
mation of quantum vortices. Astronomical observations may 
thus provide a means to diagnose different dynamical ef- 
fects, thereby allowing us to constrain (or rule out) this 
form of dark matter. In fact, BEC-CDM has previously 
been invoked to overcome the above mentioned shortcom- 
ings of standard CDM: One prime motivation for consider- 
ing BECs for CDM has been their ability t o produce galactic 
haloes with constant de nsity cores, see e.g. iGoodmanl (|2000h 
("repulsive DM") and iPeebles! (|2000l ) ("fluid DM"), since 
the corresponding profiles may then better agree with ob- 
served rotation curves of dwarf and LSB galaxies as has been 



2 This is different from the kind of self-interacting, cold dark mat- 
ter partic les referred to elsewhe r e in t he literature as SIDM, sug- 
gested bvlSper eel fc S teinhardU (2000), which we have studied in 
lAhn fc Shapiro! l l2005t) and lKoda fc Shapird ll201lh . In SIDM, the 
particle self-interaction results in 2-body elastic scattering which 
adds 'collisionality' to the otherwise collisionless CDM gas, but 
does not make a BEC or exhibit any form of macroscopic quan- 
tum coherence. The reader should henceforth avoid any confusion 
between the BEC-CDM discussed here and SIDM. 



© RAS, MNRAS 000,[TH29] 



Angular momentum in BEC-CDM haloes 3 



shown in the papers oflArbev. Lesgourges fc Salad |2003h : 
iBohmer fc Harkq |2007l ) and Bohmer, Martins, Salucci (pri- 
vate communication) . The problem of overabundance of sub- 
haloes in standard CDM is naturally resolved in some BEC 
dark matter models, since the uncertainty principle pre- 
vents the formation of gravitationally bound isolated struc- 
tures below a certain length sca le, which depends on the 
dark m atter particle mass, see e.g. lHu. Barkana fc Gruzinovl 
(2000) ("fuzzy DM"). However, the formation of large- 
scale structures involving BEC dark matter has not yet 
been studied in the same depth as for standard CDM 
and many issues are still unresolved. Some recent results 
on structure formation s tudies involving BEC-CDM ca n 
be found in the works of lUreiia- Lopez fc Guzman! 12 003) ; 
iFukuyama. Morikawa fc Tatekawal l!200ct); I Woo fc Chiuehl 
d2009l); [Marsh fc Ferreiral (|2010h ; lHarkol |201lh and 
IChavanisI (|2012| 1 3 I. 

We will see in the next section that for quantum- 
coherence to influence structure on the scale of galactic 
haloes of radius R, either the de-Broglie wavelength of the 
dark matter particle is of that same order, 



Ad 



requiring the particle mass to be of the order 

m>m H = 10" 25 (i?/100 kpc)~ 1/2 (M/10 12 M )" 1/2 eV, 

or else 



A,, 



< R 



and the halo is supported against gravity by self-interaction 
pressure, which requires both 

m > m H 

and 

g > g„ 2 ■ 10 _64 (i?/100 kpc)(M/10 12 M ) _1 eVcm 3 , 

and we choose units where c = 1. 

However, previous literature on BEC-CDM has mostly 
neglected an important aspect of halo physics, namely an- 
gular momentum. In the early phases of halo collapse, tidal 
torques caused by large-scale structure give a halo most of 
its angular momentum. Cosmological N-body simulations of 
the standard CDM universe show that haloes form with a 
net angular momentum such that the dimensionless ratio 

A- ii£|1/2 (1) 

which expresses their degree of rotational support, has typ- 
ical values in the range [0.01, . 11 with median value ~ 0.05 . 
These values can be found in iBarnes fc Efstathioul (|l987t ). 
and more recently in lAntonucci et al. (2010) , which s eem t o 
be confirmed by observations, see iHernandez et al] |2007). 



3 ISlepian fc Goodman! 1120111 ) recently suggested that astronomi- 
cal constraints can be used to rule out strongly-repulsive bosonic 
dark matter. However, their analysis does not apply to the case 
considered here. They assume that the bosons are in thermody- 
namic equilibrium in isothermal haloes at the halo virial temper- 
ature, with a condensate core surrounded by a non-condensate 
envelope. As they themselves point out, their assumption breaks 
down when 2-body collisions are not frequent enough, which is 
the case here (see also Section 5). 



In the above expression, L is the total angular momen- 
tum, E is the total energy, and M is the total mass of 
the halo. The quantity A 2 corresponds roughly to the so- 
called i-parameter, t = T/\W\, where T is the rotational 
kinetic energy and W is the gravitational potential energy. 
Tidal-torque theory can successfully account for the A- values 
found in A^-body simulations of stru cture formation, see e.g. 
iPorciani. Dekel fc Hoffman! (j2002al fbh. We shall here be in- 
terested in the case where the BEC nature of CDM affects 
small-scale structure and the internal dynamics of galactic 
haloes, while large-scale structure formation follows that of 
standard CDM to a great extent. In this paper, therefore 
we will adopt the above range of A-values for BEC-CDM 
haloes, too. The general problem of acquisition of angular 
momentum by BEC haloes is worth further consideration in 
a future work, but in the context of this work, we will content 
ourselves with the convenience of the above prescription. 

Once haloes rotate, additional effects come into play. 
Halo shapes and profiles will differ from those which result 
when angular momentum is zero, and, since self-interacting 
BEC-CDM behaves as an irrotational superfluid, these will 
differ from the case of collisionless CDM. In addition, it 
is known that laboratory BECs can develop quantum vor- 
tices when rotated wi th a sufficient angular velocity (see e.g. 
iMad ison et al] (2000)), thereby changing the density profile 
in a characteristic way. The question arises as to whether 
this may be also possible for BEC haloes and we will at- 
tempt to answer that question here. 

Se veral authors have p r eviou sly pursued related ques- 
tions. ISilverman fc Mailed (|2002l ) postulated vortices in 
galactic haloes comprised of ultra-light bosonic particles by 
comparing the critical angular velocity for vortex forma- 
tion with the rotation rate of M31. However, the formula 
quoted there without derivation applies actually to labo- 
ratory superfluids with st rong self-interaction but without 
self-gravity. Subsequently, lYu fc Morgan! (|2002h have stud- 
ied the influence of a vortex lattice on the velocity profile of 
a spherical galactic hal o composed of ul t ra-lig ht bosons. 

On the other hand. lDuffv fc Sikivid (|2008l ) have argued 
that certain fine-structure in the observed inner mass dis- 
tribution of the Milky Way can be explained only if the 
infalling dark matter particles had a net overall rotation, 
causing a 'tricusp' caustic ring in the catastrophe structure 
of dark matter in that case. For standard, non-interacting 
CDM models, however, one expects infall to be irrotational. 
According to ISikivie fc Yangl lj2009t ). axionic dark matter, 
as a BEC, may be able to form a vortex lattice with high 
enough vorticity so as to mimic a net rotational component, 
producing the above structure. As such, the authors suggest, 
Milky- Way observations may already have detected the sig- 
nature of axionic dark matter. We will come back to this 
claim in Section 5. 

This paper will extend the ana lysis which we presented 
in iRindler-Daller fc Shapiro! (|2010h . We showed there that 
vortices are favoured to form for a wide range of possible 
particle mass and self-interaction strength, by calculating for 
the first time the critical angular velocity for vortex creation 
for a simple model of BEC-CDM galactic haloes and com- 
paring the result with the angular velocity expected from 
cosmological N-bod y simulations of standard CDM haloes. 
iKain fc Lind (|2010f ) subsequently studied the formation of 
a vortex in spherical (non-rotating) haloes in the Thomas- 



© RAS, MNRAS 000,[TH29] 



4 T. Rindler-Daller, P.R. Shapiro 



Fermi regime of strong self-interaction. The major draw- 
back of their approach is the fact that it starts from haloes 
which have no angular momentum. The physical mechanism 
by which vortices shall form is thus absent from the begin- 
ning. In fact, the amount of angular momentum of a singly- 
quantized, axisymmetric vortex in the center of a halo is 
given by 



Nh, 



(2) 



and this angular momentum has to be provided by halo spin- 
ning. We will report here our study of the equilibrium of ro- 
tating, self-gravitating, virialized haloes, and the conditions 
for which vortex formation is favoured. For the latter, we 
will focus our attention on two extreme cases, one in which 
haloes have just enough angular momentum to support a 
single vortex, i.e. their angular momentum is L = Lqm, and 
the other case in which haloes have much excess angular mo- 
mentum, such that L ;§> Lqm- Haloes will be modelled in 
each case as rotating, ellipsoidal bod i es. W hile the main con- 
clusion of iRindler-Daller fc Shapiro! (|201dh will not change, 
this paper will significantly improve on the analytic mod- 
eling of BEC haloes presented there and thereby establish 
better constraints on the BEC dark matter particle mass 
and self-interaction required for vortex formation to hap- 
pen. We will employ an energy analysis of the full equations 
of motion, thereby deriving analytically the conditions un- 
der which vortex formation lowers the energy and is, hence, 
favoured in haloes with the amount of angular momentum 
expected from large-scale structure in the CDM universe. 

The impact of vortices may be profound: We can expect 
vortices to reside preferentially in the centers of dark mat- 
ter haloes. Vortices deplete the dark matter density in their 
core region once they have formed. This in turn influences 
the dark matter density profile and also the gravitational 
coupling to the baryons. If less dark matter is around than 
without vortices, for example, then baryon cooling and con- 
densation, a prerequisite for star formation, may as a con- 
sequence be reduced and delayed. It is thus interesting to 
ask for which BEC-CDM models such vortices can be ex- 
pected to form. However, we stress that this work will be a 
'dark matter only'- analysis, and the baryonic component of 
galactic haloes will not be considered. 

This paper is organized as follows: In Section 2, we re- 
view the basic equations, the self-interaction regimes of BEC 
dark matter and their consequences for halo density pro- 
files. The underlying equations we are going to use are the 
Gross-Pitaevskii equation for the wave function of the Bose- 
Einstein-condensed halo, and the Poisson equation by which 
the density is coupled to the halo gravitational potential. 
The equation of motion can be restated in terms of quantum- 
mechanical fluid equations. So, in Section 3 we will study 
approximate figures of equilibrium models for BEC haloes 
with angular momentum included: rotating haloes will be 
approximated either by irrotational Riemann-S ellipsoids or 
by Maclaurin spheroids. In Section 4, we will present the en- 
ergy analysis to determine when vortex formation will occur 
in such BEC haloes, whose net amount of angular momen- 
tum we will fix by values for the A-spin parameter, which 
are representative of those of standard CDM. In the regime 
of strong self-interaction, we will thereby establish bounds 
on the BEC dark matter particle mass and self-interaction 
strength above which vortices will be favoured. The implica- 



tions of our results for halo models comprised of BEC-CDM 
will be presented in Section 5, which contains our conclu- 
sions and discussion. Some further detailed derivations, ta- 
bles and frequently used relationships are deferred to three 
appendices. 



2 FUNDAMENTALS OF BEC DARK MATTER 
HALOES 

2.1 Basic equations 

We will describe self-gravitating BEC haloes by self- 
consistently coupling the Gross-Pitaevskii (CP) equation of 
mot ion for the complex sca l ar BE C wavefunction ip(r, t) (see 
e.g. IPitaevskii fc Stringaril ||2003h ) of the dark matter halo 
to the Poisson equation, 

ih ^ = -^ A ^ + ( m $ + 9\4>\ 2 )4>, (3) 

A$ = 4-irGm\iP\ 2 , (4) 

where |^| 2 (r, t) = n(r, t) is the number density of dark mat- 
ter particles of mass m and $(r) is the gravitational poten- 
tial of the halo. We assume that all N particles comprising a 
given halo of volume V are in the condensed state described 
by ip, such that 



N. 



(5) 



This effectively amounts to assuming that the gas is at zero 
temperature. Self-gravitating BEC matter has been con- 
sidered before as a possible candidate for making ultra- 
dense stars, called bos on star s , that can avoid collapse to 



blac k hole (see e.g. Kaurj (119681) : i Ruffini fc Bonazzolal 
l| 19691 ): IColpi. Shapiro fc Wassermanl (|l986l )). BEC haloes 
can be imagined as boson stars on galactic scales, except 
that their densities are very low. General-relativistic effects 
are thus usually considered to be negligible on those scales. 
As such, we are free here to treat the effects of gravity in the 
Newtonian limit, as is customary in the literature on BEC 
haloes. 

BEC-CDM, like standard CDM, is assumed to interact 
so weakly with other matter and radiation, once its abun- 
dance is fixed in the early universe, that we can neglect all 
other, non-gravitational couplings. However, in contrast to 
standard CDM, BEC dark matter can be self-interacting. 
The BEC self-interaction has been included in equ.(j3]l in 
the usual way in terms of an effective interaction potential 
g\i/j\ 4 /'2 with coupling constant (or self-interaction strength) 
g. The possibly complicated particle interactions are simpli- 
fied this way in the CP framework in the low-energy limit of 
a dilute gas: disregarding higher than 2-body interactions, 
the cross section for elastic scattering of indistinguishable 
particles becomes constant in the low-energy limit, 



■ 87ra„ 



(6) 



with the s-wave scattering length a a . The coupling constant 
of the effective interaction is then given by 

g = 4nK 2 — . (7) 

In this work, we will restrict our consideration to g ^ 
(a s ^ 0). Condensates having attractive particle interactions 
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g < have been shown not to support vortices, but rather 
produce bright solitons with which we will not be concerned 
here. We shall note that the above GP equation is strictly 
valid for dilute quantum gases only, which means that a s 
must be much smaller than the mean interparticle distance, 
i.e. a s < n~ 1/3 . 

The complex wave function of the condensate in the GP 
equation, a form of non-linear Schrodinger equation ((3]), can 
be decomposed into its amplitude and phase function, 

^r,t) = \i>\(r,t)e iS ^ = J^e i3 ^ (8) 

V m 

with the corresponding dark matter halo mass density 

p(r,t) = m\4,\ 2 . (9) 

Inserting (|8]) into (|3}, the GP equation decouples into two 
equations for the real functions \ip\ and S, 



2m 
and 
dt 



dS_ 

dt 



+A|t/>|-|i/>|(VS') 2 



2m 
1?~ 



(m$+3|Vr)|V>l = 0(10) 



IVl 2 -vs 

m 



= 0. 



The associated quantum-mechanical current density, 



JL(^*V^- il>Vip*) = n{v,t)— VS(r,t), 
Aim m 



(11) 



(12) 



can be expressed in terms of the bulk velocity v of the gas, 
if we write 



vs. 



(13) 



The GP equation <[3j) can hence via equ. (|10|l and (|lip be 
written as a system of quantum-mechanical hydrodynamic 
equations for the mass density p and the velocity v, easily 
recast in the form of an Euler-like equation of motion along 
with the continuity equation, 



0~v 

P-qI +P(v- V)v = 

| + V.(pv) = 0, 
where we define 

h 2 A^p 



-pVQ- P V<S>-VP S i, 



Q 



2m 2 yfp 



(14) 
(15) 

(16) 



The term Q gives rise to what is often called 'quantum pres- 
sure', an additional force on the rhs of equ. (|14|l which basi- 
cally stems from the quantum-mechanical uncertainty prin- 
ciple. The particle self-interaction, on the other hand, gives 
rise to a pressure of polytropic form 



Psi = K p p ' 



9 

2m?' 



(17) 



where the polytropic index is n — 1 (not to be confused 
with the number density n(r, t)), and where the polytropic 
constant K p depends only on the dark matter particle pa- 
rameters. According to equations (|14[) - (|17[) , BEC haloes as 
quantum gases have fluid-like properties, in contrast to the 
collisionless nature of standard CDM haloes. 

The definition of the bulk velocity in equ. (|13[) implies 
that for any smooth phase function S, the velocity flow of the 



system is irrotational, V x v = 0, which is a typical charac- 
teristic of superfluids. However, in the presence of a vortex, 
the fluid velocity experiences a singularity and the circula- 
tion around a contour enclosing the vortex is non-vanishing 
and a multiple integer of the elementary circulation, 



VS • dl 



2-nd— 
m 



(18) 



where dl is a unit tangent vector to the curve C encircling 
the vortex with winding number d. The parameter d is an 
integer in order to ensure that the wave function, whose 
amplitude vanishes along the vortex, is singly-valued. The 
genuine quantum character of the system can be inferred 
by the fact that the above fluid circulation condition has 
effectively become a 'quantization condition'. For an axial- 
symmetric vortex, S = dcj>, the velocity flow around it is 

v=A^. (19) 
m r 

A wavefunction having such a vortex is an eigenstate of the 
angular momentum with l z = dh, such that the vortex car- 
ries a total angular momentum equal to 



dNh = dLQM- 



(20) 



Lqm, as already defined in is the minimum angular 
momentum necessary to sustain a singly-quantized vortex 
and will constitute an important quantity of the analysis in 
Section 3 and 4- 

The hydrodynamic equations described above make it 
possible for us to apply familiar results for self-gravitating, 
classical fluids to derive properties of quantum-mechanical 
BEC dark matter (e.g. the classical figures of gravitational 
equilibrium, the properties of vorticity-free flow, the con- 
ditions leading to gravitational instability or vortex forma- 
tion). In this way, we can also make clear the distinct role 
played by the quantum pressure term in the Euler equation 
in contrast to the roles of the polytropic self-interaction pres- 
sure and gravity. In the future, this will also be useful when 
we consider the dynamical origin of the virialized structures 
which constitute the cosmological haloes in this model, in- 
cluding the application of numerical hydrodynamic meth- 
ods. 

In this paper, however, we shall limit our treatment to 
the equilibrium structure of the virialized haloes expected 
to emerge during the cosmological formation of galaxies and 
large-scale structure in a BEC-CDM-dominated universe. 
The dynamical formation of such structures is presumed to 
result from the gravitational instability of primordial density 
fluctuations, like those in the standard CDM model, lead- 
ing to the 'cosmic web' of filamentary structure. In order to 
follow the development of vortices, predicted here by a sta- 
bility analysis, an initial-value problem involving the fully 
time-dependent 3D evolution from linear perturbations to 
the turn-around and nonlinear collapse of individual haloes 
must be solved, including the tidal torques responsible for 
halo angular momentum. This is a challenging task and be- 
yond the scope of this paper. Here, instead, we take the first 
step of determining whether the virialized haloes in a BEC- 
CDM universe are subject to vortex formation by applying 
an energy argument to the stationary description of BEC- 
CDM haloes as rotating, self-gravitating figures of equilib- 
rium. If the energy of such a halo is lower in the presence 
of a single, central vortex than without a vortex for some 
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values of the dark matter particle mass and self-interaction 
coupling strength, for rotation rates consistent with the an- 
gular momentum expected from large-scale structure forma- 
tion, then the halo is unstable to vortex formation for those 
parameters. Before we turn our attention to rotating haloes, 
we set the stage by reviewing two important limiting cases 
which give insight into how the properties of the DM particle 
can lead to quantum-coherence on galactic scales. 



2.2 BEC dark matter halo regimes 

In this paper, we will study bound, isolated systems whose 
equilibrium structure must satisfy the Virial theorem (see 
relationship (|59|l in the next subsectionj^. BEC haloes can 
be stabilized against gravitational collapse by internal pres- 
sure or rotation. In addition, BEC haloes can be prevented 
from collapse by the quantum-mechanical uncertainty prin- 
ciple, which is expressed by the quantum-kinetic term in 
((31 . In order to gain an idea for the length scales involved, 
we may simply check the importance of terms of the time- 
independent GP equation in different limiting case sfl Ob- 
viously, the corresponding characteristic length scales will 
be determined by the respective particle or halo parameters 
which are involved in the terms under consideration. In the 
following, let R denote the size, M the total mass and p the 
mean density of a given halo. 



2.2.1 Non-interacting DM particles: quantum pressure 
versus gravity 

For quantum-coherence to be relevant on the scale of a halo 
the dark matter particle de-Broglie wavelength 

h 



X 



(21) 



should be of the order of the size of the system 

X deB < R. (22) 

Using v V i r ~ Vdrc = {GM / R) 1 / 2 in the expression for XdeB, 
the particle mass must be of the order of 

m " (GMR) 1 / 2 ~ R 2 {GpY' 2 ' (23) 

For galactic virial velocities of the order of v V ir ~ 10 — 200 
km/s, it is clear that it is the particle mass which must be 
small enough in order for the de-Broglie wavelength to be 
comparable to the extent of the halo. 

The quantum-kinetic term can stabilize a system 
against gravitational collapse, 



= Aib + m$ib 

2m 

above a length scale of 



(24) 



4 The dynamical evolution which leads to the formation of the 
virialized haloes in a BEC-CDM cosmology is beyond the scope 
of this paper, which focuses on the haloes themselves. For some 
di scussion of this d ynam i cal evolution, the reader is ref erred 
to lAlcubierre et al.l d2002h : iGuzman fc Urefia-Lopezl {2004 ) and 
IShort fc Colesl d2006t) . 

5 A Jeans instabilit y analysis in the se limiting cases has been 
recently presented m lChavani 1 tecml) . 



which 
(see 



h 



m(2Gpy/ 2 



plays 
also 



the role of a quantum 



(25) 

Jeans lengt h 
(1985); 



■Khlopov. Malomed fe Zeldovichl 
iHu. Barkana fc Gruzinovl (|2000h '). In BEC-CDM with- 
out self-interaction, Iqp is the smallest scale above which 
bound structures can form. More precisely, it can be shown 
that the density profile of a bound system like a halo, can 
on ly be determined numerically, and is given for instance 
in iMembrado. Pacheco fc Sanudg ()l989l ) . It has a flat core 
and falls off as r~ 4 at infinity. Although the total mass 
M — Nm is conserved and hence finite, the system has no 
finite size (i.e. no compact support) in this regime, and so 
one may calculate a radius which includes 99 % of the mass, 



#99 = 9.9 



h 2 



GMm 2 



(26) 



IMembrado. Pacheco fc Safiudel l|l989l ). Indeed, this radius 
which is about 5 times as large as (|25[) . is just the more 
accurate form of Iqp, which can be seen by expressing the 
density in (|25[) by using the total mass and radius, instead. 
Now, if we require XdeB to be of the order of 



Iqp i; XdeB < R, 



(27) 



the particle mass has to satisfy m < V 187r 3 mfr, where we 
have defined 



h 



R 2 {nGp) 1 / 2 



2h 



(RM) 



-1/2 



(28) 



depending on halo properties. This mass of a non-interacting 
DM particle is hence characteristic for making the dark mat- 
ter of a halo quantum-coherent on a substantial fraction of 
its size. 

There is, however, yet another meaning for this defini- 
tion, which is of importance in the context of our paper: If 
a halo of mass M and angular momentum L = Lqm = Nh 
rotates uniformly with an angular velocity SIqm, then 



Lqm = MR Q.QM 
implies 

h 



QM 



iR 2 ' 



(29) 



(30) 



On the other hand, the characteristic gravitational angular- 
frequency is usually defined as 



\J nGp. 



(31) 



It is basically the inverse free-fall time of a self-gravitating 
body, i.e. this angular velocity corresponds to a rotation, 
which supports the halo against gravitational collapse. Now, 
the dark matter particle mass assumes the form (|28[) if 
£Iq = Qqai, i.e. if the halo's angular velocity due to gravita- 
tional stability is sufficient to have each particle contribute 
an amount of h to the total angular momentum. 

This regime, where self-interaction is neglected, is de- 
noted as Regime I in Table [T] 
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2.2.2 Strongly-interacting DM particles: self -interaction 
pressure versus gravity 

In the regime where the particle self-interaction dominates, 
the so-called Thomas-Fermi (TF) regime^, it is only the self- 
interaction which balances gravity, 



= m<f>V> + g\tp\ tp. 



(32) 



A system is then gravitationally bound above a length scale 
of 



Is 



(33) 



More precisely, the GP equation of state is exactly an 
(n = l)-polytropic law in this case, and equ. (|14[) reduces to 
a Lane-Emden equation of motion. Solving (|51[) with (|52[) . 
or equivalently (|14|l in this regime, leads to the following 
general result for the mass density profile, 



in 
9 



P - 



m o 



(34) 



which can be determined, once $ and v are known. For 
spherical haloes in the static case, the density profile is given 
by 



P (r) 



s . 

Pc sin 



2-kG 



2-kG 



(35) 



with the corresponding ratio of central to mean density as 

(36) 



pf lp S = tt 2 /3 ~ 3.29 



A spherical BEC halo is thus an (n — l)-polytrope with 
radius 



Ro 



2ttG 



AirGm 2 ' 



(37) 



This result, derived in IChandrasekhail (|l939h for stellar 
bodies, has been re-derived in t he context of long-range- 
interacting atomic BEC gases in lO'Dell et al.1 ^2000) and 
in the context of d a rk matt e r, for insta nce in the works 
of iGoodmanl l|200Ch . iPeebled (|2000l 1 and iBohmer fc Harkol 
(|2007h . The Thomas-Fermi regime will be studied in more 
detail in Section 3. Since this limit corresponds to the BEC 
becoming an (n = l)-polytrope, we will refer to this case 
also as the polytropic regime. The convenient simplifications, 
which come by disregarding the complicated quantum pres- 
sure term in (|14p in this regime, have been heavily exploited 
in much of the literature on the subject. However, this comes 
with a price. According to equ. ((3"7)) . the size of the poly- 
trope is fixed by the BEC-CDM particle mass and coupling 
strength in the combination g/m 2 . While we shall continue 
to refer to this size as the halo size, it shall be understood 
that this refers either to the actual halo size or to the size of 
a virialized core region within a larger halo. The TF regime 
considered here is denoted as Regime II in Table 1. 



6 The Thomas-Fermi energy functional appeared originally in 
atomic physics as an energy depending only on the density of 
the system. Since the GP energy, too, depends on density only in 
the above limit case, its formal analogy has been given the same 
name. 



2.2.3 Validity of the Thomas-Fermi regime: quantum 
pressure versus self-interaction 

To determine whether a halo is in Regime I or 77, we must 
compare the quantum pressure and self-interaction pressure 
terms to each other. When the length scale for density vari- 
ation in the BEC fluid is such that the two terms are of 
equal magnitude, that length scale is sometimes called the 
'healing length'. The reason for this name is because it 
is the characteristic distance over which the wavefunction 
approaches the background value of a smooth density dis- 
tribution, subjected to a localized perturbation. Balancing 
quantum-kinetic term and self-interaction, 

= -|-A^ + 5l^| 2 V', (38) 

results thus in the following expression for the healing 
length, 

h 



(39) 



We will see below that if the system's characteristic size is 
much larger than this, R 3> £, tharQ the quantum pres- 
sure can be neglected compared to the self-interaction term 
and we are in the Thomas- Fermi regime (Regime II). For 
a given mean halo density p, we see that £ grows with de- 
creasing coupling g, approaching the opposite regime of non- 
interacting particles. On the other hand, if g is fixed, the 
healing length is smaller for haloes of higher mean density. 
Let us define the characteristic coupling strength for which 
R = f, therefore, as gn, given by 



h 2 



9H 



2pR? 



2 2 R 

r h m' 



(40) 



where gu is determined by the mean density and radius 
(or mass and mean radius) of the halo. In later sections, it 
will be convenient to describe dimensionless BEC particle 
parameters, and so already here we define the quantities 



V 



and x = . / — = — . 



(41) 



Thus, the BEC self-interaction regime can also be deter- 
mined by the ratio of g /gn ■ 

Let us investigate this in more detail: In the TF-regime, 
the quantum-mechanical pressure in (| 14p is neglected. In or- 
der to see for which part of the BEC parameter space this 
regime can be expected to be a good approximation, we 
estimate the respective forces associated with the quantum- 
pressure term and the self-interaction term in (|14[1 . respec- 
tively. Both depend on the halo density profile's varying on a 
length scale of order R. Using (|16|) . the force due to quantum 
pressure is estimated by 



(42) 



2m 2 R 3 ' 

while the force due to the self-interaction is of order 

VPsi^-tMt;- (43) 



2m? R 

The ratio of the two becomes 

7 Since the diluteness condition implies n~ 1//3 <C £, one should 
consider the limit R/£ — ¥ oo, rather than £ — > 0. 
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h 2 



-pVQI 

VPs/j ~ g P R? 



,9h 



(44) 



by using the definition of gn- We see that for the above ratio 
to be much smaller than one, we require that 



9H 



>2, 



(45) 



which is thus the condition for the TF regime to be valid. 
If we go back to the characteristic size of a spherical halo in 
(|37[) . we can furthermore see that, if we set v = v vir in (121[) , 
we can rewrite 



^0 — — ( ^deB, 

12 \gHj 



(46) 



from which follows R /X deB (#/ff£r) 1/2 = Ro/£ > 1- 

Therefore, in contrast to Regime I, XdeB is much smaller 
than the size of the halo. However, using (|41l) . we see that 
XdeB — 4.3£ in this regime, so the de-Broglie wavelength is 
still larger than the healing length. In Section 3, we will ap- 
ply the Virial theorem to ellipsoidal BEC haloes in the TF 
regime, and we will see that this enforces the particle pa- 
rameters m and g to be related according to (|81[) - (|82[ ) an d 
(|103[) - (|1U4[1 . respectively. In the spherical case, the dimen- 
sionless quantities (|4ip are then related as 



(47) 



Therefore, it follows that in the TF regime we also have 
m/rriH 3> 1. In other words, XdeB can be much smaller than 
R and hence the particle mass m can be much larger than 
ran, but only if the self-interaction is high enough. 

The above characteristic units for the DM particle mass 
and coupling strength can be written in fiducial values as 



m H = 1.066-10" 
and 

g H = 2.252 ■ 10" 



R 



100 kpc 



-1/2 



M 



10 12 M Q 



-1/2 



eV, (48) 



R 



100 kpc J \ 10 12 M© 



M 



eV cm 3 , (49) 



in units where c = 1 and [R] = kpc, [M] = Mq. For haloes 
with a size between the Milky Way (M = W 12 M Q ,R = 100 
kpc) and a typical dwarf galaxy (M = 1O 1O M ,/? = 10 
kpc), we see that the above quantities are in the range of 



10" 



eV, g H 



2 . 10 - 64.-63 e y cm 3^ 



(50) 



The two BEC dark matter regimes are schematically 
described in Table [TJ where also some common names are 
listed. The latter, however, have appeared in more recent 
papers, even though both regimes have been studied already 
in earlier literature on self-gravitating bosonic matter. 

To re-iterate, both regimes revisited in this section ex- 
hibit non-classical behaviour: in Regime I, VQ in equ. ()14[) 
causes the Heisenberg uncertainty principle to act on galac- 
tic scales. In Regime II, on the other hand, VPsi results in 
(n — l)-polytropic solutions. However, the physical nature 
of this pressure and the so-caused "internal energy" stems 
from the repulsive quantum-mechanical 2-body scattering 
at vanishing temperature. 



Table 1. BEC dark matter regimes for a halo of (mean) radius 
R and mass M 



I 


II 


quantum 2> self-interaction 


self-interaction 2> quantum 


pressure 


pressure 


£-> R, g -> 0, 


g->g H , 


m<m H = f(R, M) 


m ^> ran, g/m 2 = f(R) 


R > XdeB > iQP 


R = Ro — Isi 3> XdeB 


— > grav. bound, 


— ¥ grav. bound, 


halo core ~ XdeB 


halo core ~ Rq 


no vortices 


vortices favoured 


"fuzzy dark matter" 


"repulsive/fluid dark matter" 


(free scalar-field DM) 


(scalar-field DM with strong, 




positive self-interaction) 



2.3 Stationary systems and virial equilibrium 

Since we are going to study the energetic stability of vortices 
in this paper, we will restrict ourselves to some simplified 
models in the forthcoming analysis: First, we will consider 
stationary systems and their corresponding energy, i.e. we 
will restrict to time-independent systems. Our second major 
assumption will be that haloes rotate with a constant angu- 
lar velocity fl about their rotation axis. This will allow us 
later to move into the frame co-rotating with the halo and 
to study stable structures in that frame. 

Stationary states are described by wavefunctions of the 
form ip{r,t) = ip s (r)e~ Hlt ^ h , where \i, the GP chemical po- 
tential, is fixed by the conservation of particle number. For 
these states, the mass density p — m\ip s \ 2 and, hence, the 
gravitational potential $ are time-independent, while the 
wavefunction evolves harmonically in time. Inserting this ij) 
into ((3| results in the time-independent GP equation with 
eigenvalues /i (see also Appendix C), 



h 

-—A + g\ip s \ 2 + m$ ) ips = iiips 
2m 



(51) 



The time-independent part i/'s( r ) itself can be decomposed 
as 



Mr) = \4>.\(i 



(52) 



with amplitude and phase both depending on position only. 
We will omit the subscript 's' in the forthcoming analysis. 
Systems obeying (|5ip can be studied via the corresponding 
GP energy functional, which is given by 



em = 



w 4 



j3 

a r. 



(53) 



On the other hand, equ. (|5ip could have been also obtained 
by variation of (|53p with respect to (or its complex con- 
jugate ip*), under the constraint ([5]), 



5S\ 



|V| 2 d 3 r = 0, 



(54) 



5i/j* 8i/j* 

with fi playing the role of the Lagrange multiplier. Inserting 
((52)) into (l53l> . the total energy can be written in a very 
instructive way as 

E = K + W + Usi, (55) 
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where the total kinetic energy term K is given by 



(56) 



with Kq accounting for the quantum-kinetic energy and T 
for the bulk kinetic energy (rotational, internal motion, ect.) 
of the body. Kq has no classic counterpart, and is absent in 
the cla ssical figures of equilibrium s t udied in lChandrasekharl 
or in lLai. Rasio fc Shapiro! (|l993l ). Also, Kq is ne- 
glected whenever the Thomas- Fermi regime is employed. The 
other terms in (|55p are simply the gravitational potential en- 
ergy 



W= I ^d 3 r 



and the internal energy 



Usi 



2m? 



2 j'i 

pdr, 



(57) 



(58) 



which is determined by the particle interactions. We have 
defined the latter essentially as Usi = J PsidV with Psi 
the pressure due to self- interaction (JTTJ) , which, in general, 
is not the total pressure given the presence of Kq. While 
we call Usi the internal energy, it shall be kept in mind 
that its origin is not due to thermal processes but solely due 
to the repulsive 2-body elastic scattering. The above energy 
contributions are those which enter the scalar Virial theorem 
of a rotating, isolated BEC halo under self-gravity, 

2K + W + Wsi = 0. (59) 



3 BEC HALOES WITH ANGULAR 
MOMENTUM 

In the previous section, we have reviewed halo density pro- 
files in the absence of rotation and angular momentum. 
However, we expect and assume galactic haloes to have 
undergone tidal torquing in the early phases of their col- 
lapse, and as such, we will model them as rotating, non- 
spherical bodies. In case of a (non-rotating) sphere (fl = 0), 
the BEC wavefunction were real and positive. However, for 
non-vanishing angular velocities, fl 7^ 0, superfluid cur- 
rents arise. As a result, the wave function is complex and 
therefore has a non-trivial phase or velocity flow, respec- 
tively, even without vortices (for laborat ory BEC exam ples 
see for instance the systems studi ed in iFetterl (| 19741 ) and 
iRecati. Zambelli fc Stringaril (|200ll n. This means, in partic- 
ular, that such a system has a non-vanishing, bulk angular 
momentum prior to vortex formation, which we will want 
to associate with that provided by tidal torquing, and char- 
acterized by the A-spin parameter defined in equ.JT]). 

So, in modelling dark matter haloes as BEC wave func- 
tions, ip — \4>\ e,,S j their velocity information is contained in 
the (real) phase function S(t, r). As long as no defect struc- 
tures like vortices appear in the flow, the condensate phase 
is a smooth function in which case we will use the notation 
So in the rest of the paper. The fluid velocity is generally 
given by v = h\7So/m, see eau. (|13[) . while in a frame rotat- 
ing rigidly with a constant angular velocity fl, it is given by 



v' = v — fl x r. In Section 4, we will study the energetic 
stability of vortices, and it will turn out to be convenient 
to move into a frame which rotates with the figure, i.e. at 
an angular velocity fl. It is this co-rotating frame we shall 
consider, assuming that the haloes rotate about the z-axis, 
such that fl = (0,0, fl). 

As we have pointed out in Section 2, the (vortex-free) 
density profile of a BEC halo has a flat core, in contrast 
to the r~ -central cusps of standard CDM. However, that 
equ. Q - with or without rotation - favours a flat core over 
a cusp can be also seen as follows: The equilibrium matter 
distribution of the BEC halo in a frame rotating rigidly with 
constant angular velocity f2 is given by 

+ (m$' + g\ip'\ 2 - ft')IV>'l ~ ?#'|VS' • (fi x r) = 0, (60) 

(see (|10p and (|51[) ) where primes will in the following de- 
note quantities as measured in the rotating frame, and 
/j,' = const. Inserting a cuspy test function, | 2 = (r')~ a /m 
with constant exponent a > 0, results in 



m $' + i.( r ')— + |-(VS' 

= li > + £.2(« + 1 ) (/)- 

P 2m 2 \ 2 J v ' 



hVS' • (O x r) 



(61) 



Calculating the associated gravitational potential $' to the 
test function by requiring that V'$' = at r' = (i.e. 
no net gravitational force at the centre) requires a < 1. 
However, if a > and since p! = const., one can easily 
show that the resulting velocity flows v' would diverge in the 
centre, regardless of the absence of a vortex, which would 
contradict the irrotationality constraint. In fact, in the next 
section, we will study figures of revolution as halo models, 
whose velocity fields prior to vortex formation are smooth 
in either frame of reference. 



3.1 BEC haloes as Maclaurin spheroids 

3.1.1 General properties 

For the sake of computational simplicity and in order to 
gain analytical insight, we start our analysis by consider- 
ing a simple halo model, the densities and potentials of 
which are given by homogeneous Maclaurin spheroids, rotat- 
ing uniformly with angular velocity fl = (0, 0, fl). Although 
the BEC fluid is compressible, as described in Section 2, 
the spherical polytrope solution shows that the density only 
varies by a factor of about three from the centre to the 
surface, see (|36[) . Hence, the assumption of a uniform den- 
sity as a first approximation is more justified for BEC than 
for standard CDM haloes. However, homogeneous, i.e. in- 
compressible Maclaurin spheroids are only approximate so- 
lutions of the CP equation (|51[) . In order for our model to 
be in global virial equilibrium, cqu. (|59p , in accordance with 
the equations of motion, we will have to find a constraint on 
the underlying DM particle parameters, as will be shown in 
the next subsection. 

We denote the mass density as po, and the semi-axes 
(ai, 012,013), lying along (x,y,z), are such that a\ = a-z = 
a > as = c. The total volume of a spheroid is V = 47ra 2 c/3, 
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and we will make use of its mean radius, denned as R = 
(a 2 c) 1 ^ 3 . The gravitational potential inside the spheroidal 
body is given by 



$o(r,2) = TrGpox 



{A x {e)r 2 + A z (e)z 2 ) - 2^1 -e 2 a 



2 arcsin(e) 



(62) 



in cylindrical coordinat e s (r, z ) (see IChandrasekharl (|l969l ) 
or iBinnev fe Tremaind (| 1987h for definitions and formu- 
las, but note the difference in the sign convention of $o 
between the two re f erence s. We choose the convention in 
IBinnev fc Tremaind i|l987h ). The functions Ai and A3 de- 
pend on the eccentricity e = y/1 — (c/a) 2 of the figure. Their 
general definition can be found in Appendix A. Here, they 
reduce to functions in closed-form, 



Ai(e) = A 2 (e) = 



■ arcsin(e) 



Me) 



2^r 



(63) 



(64) 



and both are positive for eg [0,1). A family of spheroids is 
characterized by the following relationship between angular 
velocity and eccent r icity, the Maclaurin formula (see e.g. 
IBinnev fc Tremaind (| 19871 ) ) 



V2(A 1 (e) - (1 - e 2 )A 3 (e)) 1/2 



(65) 



with in (|31|) . Relation (I65|l parametrizes the rotation via 
e, and e = for Q = 0. The figure is stab le only for values of 
e below 0.9529 (see lChandrasekharl l|l969h ). For our analysis, 
we will also need the gravitational potential and rotational 
kinetic energy of the homogeneous Maclaurin spheroid, 



3 GM 



5 R 



W = 
and 



4-7T 

T= — (1 
15 v 



(1 



2sl/6 



arcsin(e) 



-1/3 



p n 2 R 5 , 



respectively, along with its angular momentum 
L = |L| = ^(l-e^-^TMR 2 



(66) 



(67) 



(68) 



(see a gain IChandrasekharl (Il969l ) and IBinnev fc Tremaind 
|l987h for reference). The Maclaurin spheroid experiences 
rigid rotation only, therefore its velocity in the two frames 
is simply given by 

v = n x r = s/2Sl G (Ai - (1 - e 2 )A 3 ) 1/2 {-y,x,0), (69) 

and v' = 0. An illustrative plot of a velocity field in the rest 
frame can be found in Fig.l. 



3.1.2 Comparison to CDM X-spvn parameter and virial 
constraint 

In this subsection, we aim to express the spin parameter A 
in equ.{T]) as a function of the halo's eccentricity only, i.e. 
its shape. To this end, we first rewrite 



A = 



L\W\ 



1/2 



1/2 



GM 5 / 2 

According to equations (|66[) and 



(70) 



L\W\ 



1/2 



GM 5 /2 



:M(1 



with the t-parameter t 
support, given by 



\W\ 6 arcsine 



GMV2 5% /5 e 
= T/|W|, a measure of rotational 



3V1 



2e arcsin(e) 



(72) 



This relation also applies to compressible Maclauri n 
spheroids as has been shown in lLai. Rasio fc Shapiro! l|l993h . 
The Virial theorem in (|59|l allows us to write 



1/2 



1/2 



1 + 



1/2 



(73) 



where we used the fact that, for a homogeneous body, Kq = 
0. Using now U31 — PqV with constant pressure Po and 
V = 4tvR 3 /3 according to the definition of the spheroid's 
mean radius R, we get for the ratio 

2Usi 2P V _ 2 5 

2m? 



-M-e 2 y 1/B ^RGV. (74) 



T \W\t t3arcsin(e) 

It remains to establish a connection between the mean ra- 
dius R of the halo and the dark matter particle parameters. 
The 'classical', homogeneous Maclaurin spheroid, which for- 
mally is an (n = 0)-polytrope, has a pressure profile and its 
internal energy vanishes. A BEC system, on the other hand, 
possesses a formal "internal energy" due to the particle self- 
interaction (158f) . and its equation of state is such that the 
pressure is constant in case of constant density. This appar- 
ent inconsistency is due to the fact that we model the BEC 
halo as a homogeneous system. In order for the BEC halo to 
fulfil the virial theorem globally, we thus set the BEC pres- 
sure equal to the average pressure of the classical spheroid, 
i.e. the DM particle parameters must be such that the re- 
sulting halo size corresponds to the one given by the classical 
spheroid. Using the pressure distribution of a homogeneous, 
ellipsoidal equilibrium configuration, 



P(r) = iiGpl [Ami ~ Aix 2 ~ A 2 y 2 



A 3 Z + —V , 



(75) 



we calculate Usi = J P(r)d 3 r, which, for the Maclaurin 
spheroid with A\ = A2, a\ — a 2 = a, a 3 = c, results in 

8 



P(r)d 3 r = ^Aa(e)(l - e 2 ) 2/ \ 2 Gp 2 R 5 



(76) 



Setting this expression equal to Usi = gpo/(2m 2 ) J d 3 r ac- 
cording to Q58| and solving for the mean radius R, we get 



R 



15 



1/2 (KvV' 2 
3 J 4 3 (e)(l-e 2 )2/3 y / \2ttG ) 

with K p in (JT7J). In the spherical limit, e 
to 

1/2 

Ro 



(77) 



0, this reduces 



15 (Ky\ 
2 \2-kGJ 



(78) 



which is a factor of about 1.15 smaller than the radius of 
a spherical (n = l)-polytropic halo given in (|37p . Equ. (|77[l 
connects the particle parameters (m, g) via K p to the mean 
radius R such that virialized haloes fulfill a certain con- 
straint as can be seen as follows: using the above definitions 
for K p and (|40j-(|4T}, we rewrite 
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Figure 1. Illustrative velocity field of a Maclaurin spheroid in the rest frame having A = 0.05. 



2K„ 



or y - 



9B 



2K v m 2 H ' 



(79) 



(80) 



9_ = Jl 

9h m 2 H (g H /m 2 H ) 

and 

2 

Inserting now the radius ()77|) . we arrive at the following 
constraint, implied by imposing virial equilibrium, 

1/2 

x. (81) 



2/0*0 = h 



v 8A 3 (e)(l-e 2 ) 2 /3, 
The corresponding dimensional relationship is accordingly 

1/2 



MS) 



8A 3 (e)(l-e 2 ) 2 / £ 



Vg, 



(82) 



i.e. virialized haloes as rotating, homogeneous Maclaurin 
spheroids lie on a line in (log m, log <7)-space, whose slope is 
completely determined by the eccentricity e and the mean 
radius R of the halo. We will make use of this result in Sec- 
tion 4- 

For now, we continue the calculation of A = A(e). Inserting 
((77)) into CO}, we get 



1/2 



1/2 



1 + 



eAz{e){l-e 2 Y' 2 



1/2 



(83) 



arcsin(e) 

In connection with (171|) , the A-spin parameter as a function 
of the eccentricity e of the halo is finally given by 



A = 



t 1 + 



e A 3 (c)(l 



,2x1/2 



1/2 



(84) 



5\/5 e V * arcsin(e) 

with t = t(e) in (|72l) . In our case, however, we fix A and 
solve for e. In what will follow in Section 4, we will take three 
representative values, A = (0.01, 0.05, 0.1), which correspond 
to eccentricities e = (0.062, 0.302, 0.550). A table of spheroid 
parameters as a function of A can be found in Appendix A. 
For given values of the halo mass density and eccentricity, 
the angular velocity is fixed. For example, for a mean dark 
matter density of 1.6 • 10~ 26 g/cm 3 for the Milky- Way and 
A = 0.05, we have Q ~ 1.3 ■ 10" 17 rad/s, so fi/n G ~ 0.22. 

The halo model presented in this subsection has the ad- 
vantage of computational simplicity, however, the approxi- 
mation of uniform rotation means its velocity flow is not 



strictly irrotational in the rest frame as required in the ab- 
sence of vortices. The approximation of uniform density, 
moreover, does not capture the effect of compressibility, 
which results from the dependence of self-interaction pres- 
sure on density in equ. (ll7p . Therefore, we will consider a 
further model, which improves upon those approximations, 
in the next subsection. 



3.2 BEC haloes as irrotational Riemann-S 
ellipsoids 

3.2.1 General properties 

In this subsection, we relieve some of the simplifying as- 
sumptions from above by allowing the density to vary in 
space and by imposing strict irrotationality in the rest frame 
prior to vortex formation. This will allow us to take account 
of the compressibility of the BEC fluid, which corresponds 
in the Thomas-Fermi regime to an (n = l)-polytrope (see 
|T7J). It is possible to generalize the Maclaurin spheroids 
of Section 3.1 to account for this compressibility in an ap - 
proximate way, as shown by lLai. Rasio fc Shapiro! l|l993h . 
However, this would not allow us to impose the additional 
constraint of irrotationality. In general, a rotating ellipsoidal 
halo cannot be both axisymmetric and irrotational if it is 
non-singular at the origin as can be seen as follows: Suppose 
without lack of generality that the velocity field of the halo 
has the form 



v = Sl(r, z) x r = v$(r, z)9 — rQ(r, z)t 



(85) 



allowing for differential rotation in r and z. The correspond- 
ing vorticity is given by 



V x v 



oz r or 



(86) 



Requiring irrotationality V x v = is now equivalent to 
requiring r 2 f2(r, z) — constant in r and Q = constant in z. 
This means that either vg = oo at r = or, if at r — 
we impose vg = 0, then it follows that £l(r, z) — 0. The 
assumption of axisymmetry along with the constraint of ir- 



© RAS, MNRAS 000,[TH29] 



12 T. Rindler-Daller, P.R. Shapiro 



rotationality lead to diverging velocity profiles in the centre 
of the halo or to a trivial solutioijf]. 

It is thus necessary to consider non-axisymmetric ob- 
jects if BEC dark matter forces us to take into account 
the irrotationality of its velocity field along with a bulk an- 
gular momentum prior to vortex formation. There exists 
indeed a 'classical' figure of rotation which can serve this 
purpose, namely the irrotational Riemann-S ellipsoid. The 
family of Riemann-S ellipsoids with semi-axes (01,02,03) 
along (x,y,z) describes uniformly rotating bodies, as in 
the Maclaurin case, but with internal velocity fields super- 
posed, which combine with the uniform rotation to yield a 
net flow with V x v = 0, see Fig. 2. Exact solutions exist 
for Riemann-S ellipsoids only for the case of uniform den- 
sity, a limitation they share with the Maclaurin spheroids 
( Chandrasckh ar] JlQQQl)^ . Fo rtunately, as we shall see below, 
lLai. Rasio fc Shapiro! (1 19931 ). abbreviated LRS93 in the fol- 
lowing, developed their 'ellipsoidal approximation' also for 
Riemann-S ellipsoids, so we shall be able to consider the 
compressible case here as weljf]. 

Again, we assume the uniform rotation to be about the 
z-axis with angular velocity fi. The internal velocity field 
having angular velocity A is required to have its uniform 
vorticity parall el to Q and to leave the ellipsoidal figure 
unchanged (see IChandrasekharl (Il969l ) and Appendix A for 
more details). The mean radius of the ellipsoid is given by 



i?= (aia 2 a 3 ) 1/3 = a 1 (l-e 1 ) 1/0 (l 
with eccentricities 



2x1/6, 



(87) 



ei = \J\ - (a 2 /ai) 2 and e 2 = s/l - (a 3 /ai) 2 . 

As in IChandrasekharl l| 19691 ) and LRS93, we write the veloc- 
ity field in the frame co-rotating with the figure as 

v' = Ciyx + C 2 xy (89) 

with 



Ci 



a? 



Ol 



2 
«2 



of + a% 0,2 af + 02 



c 



(12 



A. (90) 



(91) 



The vorticity in the rotating frame is defined as 
2 1 2 

C = (V x v ) z = A = nf R , 

a±a 2 

and the vorticity in the rest frame is given by 

C = (V x v) z = C' + 2fl = (f R + 2)a (92) 

There is a whole class of Riemann-S ellipsoids, which differ 
in their value of fa. However, only for those with /r = —2 
does the vorticity vanish in the rest frame, the so-called ir- 
rotational Riemann-S ellipsoids, and we will only consider 
those in the forthcoming analysis. It turns out that these 
ellipsoi ds fulfill ai ^ 03 ^ a 2 , i.e. they are all prolat«[3 bod- 
ies fsee IChandrasekharl (jl969)). Furthermore, we can write 



8 F or fully general-relat i vistic , rotating boson stars. iRvanl (119971 ) 
and lYoshida fc E riguchi ( 19970 assumed axisymmetry but in their 
case the velocity was singular at the centre. 

9 As shown in LRS93, the approximate solutions for the com- 
pressible case, which result from their ellipsoidal approximation, 
agree well with the true equilibria. 

10 As a curious note, we remark that, interestingly, recent studies 
on standard CDM halo fo rmation suggest that most haloes are 
of prolate shape, see e.g. iPIionis. Basilakos fc Ragone-Figueroal 



the relationship between f2 and A for fa — —2 according to 
LRS93 as 

n = 1 ( ai 1 a2 V - 2 r el * (93) 

The connection between angular velocity and axis ratios or 
eccentricities, respectively, can again be stated in a concise 
form as in 1)65[1. We take the formula from LRS93 for fjt = 

-2, 



n ■ 



n _ (2B 12 ^ 1/2 



2_Bis 



V 9" 



1/2 



1 + 



4a 2 a| 



4(l- ei ) ^ 1/2 
(2 - e\ y 



-1/2 



(94) 



with B12 — A2 — a'{(A\ — A 2 )/(a 2 — a 2 ) and the constant 
q n , depending on the polytropic index n, can be found in 
(I AI 1 1) in Appendix A. In what follows, we will again take 
advantage of formulae which have been derived in LRS93. 
These are again quantities entering the Virial theorem: the 
gravitational potential energy 



H--- . 3 GA f f(e 1 ,e2) 



5-n R 



(95) 



with 



+ 



/( ei , e2 ) = i( J 4 1 (l- e 2 )- 1 / 3 (l-e 2 )- 1 '' 3 

+A 2 (1 ~ e 2 ) 2/3 (l - e ^)" 1/3 + A 3 (l - e 2 )- 1/3 (l - e^) 2/3 ) (96) 

and Ai, A2, A3 as in Appendix A. Furthermore, we will need 
the total angular momentum 

L = ^-M (n(o? + a 2 2 ) - 2aia 2 A) z (97) 
and the rotational kinetic energy 

T = (01 ~a 2 ) 2 (tt + A) 2 + ^M( ai + a 2 ) 2 (fi- A) 2 , (98) 

with n n being another constant depending on the polytropic 
index n (see IA10|) . Neglecting the quantum pressure by re- 
stricting to the Thomas-Fermi regime will amount to setting 
n = 1 in the above formulae. 

The irrotational Riemann-S ellipsoid experiences an in- 
ternal motion on top of the rigid rotation. Using (|89|) - (|94[) . 
its velocity field in the rotating frame can be derived as 

v' = ^L_n{y,-{l-e\)x,Q) = 



2-e 2 
2H G 



/2B 12 ^ 1/2 



(8(l-ei) + ety 1/2 (y,-(l~ei)x,0),(99) 



while the velocity in the rest frame is accordingly 



e? 



v = v' + ftxr = n_ri_(y,a5,0) = 
2 e 1 



= n c 



(^r) (1 + 8(1 ~ e " )/e ' rl/2(y ' x ' 0) - (100) 



l|2006h and lGottlober fc Yepesl 1120071 1. The reason there, however, 
is to be found in details of the tidal interactions, and not the 
constraint of irrotationality. 
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Both expressions depend on the polytropic index n via q n . 
Illustrative plots can be found in Fig. 2. Note that the limit 
ei = 0, e 2 = e does not reduce to the Maclaurin spheroid 
case. Instead, the irrotational Riemann-S sequence bifur- 
cates from the non-rotating sphere. However, the sphere con- 
sidered as the first member of the irrotational sequence is 
viewed from a frame rotating with angular velocity SI, which 
can be explicitly seen by setting e\ = in the above expres- 
sions for the velocity, resulting in v = and v' = — SI x r 
(see also LRS93 and their equ.(5.5) and (5.6), which reduce 
to L = and T = 0, hence v = in that case). 



3.2.2 Comparison to CDM X-spin parameter and virial 
constraint 

In what will follow in the next section, we will also use 
the Riemann-S ellipsoid as a model for rotating haloes in 
the Thomas-Fermi regime. The condition of Virial equilib- 
rium relates the values of y = raj ran and x = \fgJg~H 
then to (n = l)-polytropes. For the general case of rotat- 



ing haloes, we use again (|79[) and (|80[1 along with formula 
(3.25) of LRS93, which relates the mean radius R of the 
ellipsoid to the radius Rq of the equilibrium spherical poly- 
trope, equ. fl37}. For (n = that formula becomes 



R = Ro[f(ei, e 2 )(l - 2t)]- 1/2 = R g(ei, e 2 )- 1/2 , 



(101) 



where the t-parameter and g = /(l — 2t) depend on the 
eccentricities of the rotating figure only, 

f(eiie2) = ^n 2 (i- e ?)- 1 / 3 (i-el)- 1 /3 x 

5 



r(2 



2\ 

ei) 



4(1 -el) 



+ 



(2 



\f(e 1 ,e 2 )\- 1 (102) 



with SI from (|94)l . /(ei, e 2 ) in (|96) and ki in (|A10|) . Inserting 
(|101[) into ((80)) and the result back into (f79|) yields the Virial 
constraint 



/ \ ^" / \ — 1/2 

= ^|ff( e l' e 2) 

with corresponding dimensional form 

t \ I" / n-1/2 "Iff /- 

V8 -Jgn 



(103) 



(104) 



depending again only on the eccentricities (ei, e 2 ) and mean 
radius R of the halo. Equ. (|103[) an d (|104() are the equiva- 
lent relations to (|8ip and ((82}, constraining the BEC dark 
matter particle parameters such that rotating haloes as irro- 
tational Riemann-S ellipsoids in the Thomas-Fermi regime 
fulfill virial equilibrium. Their formal equality does not come 
as a surprise, since the assumption of homogeneity in Sec- 
tion 3.1.2 results in the same neglect of Kq in (|59|l as in the 
Thomas-Fermi regime (albeit the latter does not neglect the 
spatial variation of the density, even though its Laplacian is 
likewise disregarded, as in the homogeneous case). 

As in Section 3.1.2, we shall now also derive a relation- 
ship between the spin parameter A and the eccentricities of 

11 For (n = 0), their formula (3.25) reduces to R = Ro. However, 
it would not have been correct to use this result for the BEC 
haloes we were considering in Section 3.1. due to our matching of 
the BEC pressure with the average pressure of an homogeneous 
spheroid, and so we needed to derive 1771 1. 



the ellipsoid (ei,e 2 ). Using l|95|>, ((97)) and J98J) for n = 1, 
the corresponding ratio 2Usi/T in (|73() is here given by 

2U S i 2U SI 



T 



\W\t 



= ifS^ 61 ^" 1 = S^-r 1/2 /Ke2)-\(105) 

and we have also used (|101|l and (|37|) . pf denotes the central 
density of the ellipsoid. After some more algebra, we finally 
arrive at 

\ 3Kl fl „2n-1/3/, 2n-1/3/ 2\ /, 4(1 — el) 



x |/(ei,e 2 )| 1/2 n [t + 



i/2o / , 2 /( e i> e 2) 



-3/2 \ V2 



3(1 - 2t)V2 



(106) 



with Cl in (|94l) and t in (|102() . We use this equation, along 
with (|A9[) . in order to solve for the eccentricities (ei,e 2 ) at 
given A. Although highly nonlinear, involving many trigono- 
metric and elliptic functions, equ. (|106|) can be solved, in 
principle, as straightforwardly as the much simpler rela- 
tionship (JHIJ. The same values of A = (0.01,0.05,0.1) as 
in the previous section correspond now to eccentricities 
of (ei,e 2 ) = (0.707,0.573), (ei,e 2 ) = (0.881,0.797) and 
(ei,e 2 ) = (0.934,0.887), respectively. Note that even for 
small A- values, the eccentricities are quite large, and sur- 
pass those in the case of the homogeneous spheroid. This is 
in part due to the fact that compressible bodies allow for 
higher eccentricities at a given angular momentum than in- 
compressible ones do. A table of ellipsoid parameters as a 
function of A can be found in Appendix A. 



4 INSTABILITY OF ROTATING BEC HALOES 
TO VORTEX FORMATION 

4.1 Energy argument 

We shall use an energy argument in order to derive the crit- 
ical angular velocity and energy for vortex creation in a ro- 
tating, self-gravitating BEC halo. Configurations rotating 
at an angular velocity Q in the rest frame will be station- 
ary solutions in the co-rotating frame. Therefore, we seek 
for vortex solutions which are (energetically) stable in the 
rotating frame. To this aim, we will compare the total en- 
ergy with and without a vortex as measured in this frame. 
The equation of motion in the rotating frame is given by 
equ.Q with an additional operator V ro t = — ■ L' on the 
right-hand-side under the brackets with angular momentum 
L' = —ihr' x V', and it is also understood that the primed 
variables and quantities are then those in the rotating frame. 
The GP energy functional is then 



-ihii>'*v'ii>' ■ (n x 



n IV7' /'|2 

2m 



+f« 2 +§iv>r- 



(107) 



In what follows, we shall use (|107[) to determine at which 
energy or angular velocity respectively, the presence of a 
vortex starts to be energetically favoured, and to derive the 
vortex energy as a function of the dark matter particle pa- 
rameters. 
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Figure 2. Illustrative velocity fields of an irrotational Riemann-S ellipsoid with n 
and in the co-rotating frame (right-hand-plot) . 



1 and A = 0.05 in the rest frame (left-hand-plot) 



We will refer to a system without vortices (thus prior to vor- 
tex formation) as the unperturbed system with ip' = f'e lS ° , 
and the unperturbed (vortex- free) halo density IV'ol 2 = 
|/' | 2 = |/| 2 . The associated halo mass density and gravi- 
tational potential will be denoted as po = ti|/| 2 and <£>o, 
respectively. The corresponding phase S' has no singularity 
in that case. These unperturbed haloes will be modelled by 
the ellipsoidal figures described in the previous section. How- 
ever, before we turn our attention to the energy analysis, we 
highlight a necessary condition for vortex existence. 



4.2 Necessary minimum condition for vortex 
formation L ^ Lqm 

Applying a rotation with small enough, finite angular veloc- 
ity to a perfectly spherical-symmetric BEC does not elevate 
its angular momentum above zero. It is only above a critical 
value Q c , when a (singly-quantized) vortex starts to form, 
that the total angular momentum is given by the amount 
necessary to sustain this vortex, Lqm = Nh. However, non- 
spherical bodies carry already a bulk angular momentum 
prior to vortex formation which is responsible for their de- 
formation, i.e. L grows with even for Q < Q, c (a very 
ill ustrative plot of such an L — Q, relationship can be found 
in iFetterl l| 19741 )). We can easily derive a relationship be- 
tween the angular momentum L of our ellipsoidal haloes and 
the minimum angular momentum Lqm necessary to sustain 
one vortex. In dividing equ. (|68[) by Lqm we have for the 
Maclaurin spheroid 

= (1-e 1 MR 2 Q = -(1 - e 2 ) 1/3 

L QM hM5 y ' 5 V ' m H ttG 

and by using (I65|l . 

T^ = ^T ^n 2 ^V3 VM*) ~ ~ e2 )^(e)- (108) 
Lqm mn 5(1 — e^) 1 '* 

The equivalent relationship for the irrotational Riemann-S 
ellipsoid using ((97]) and (J93J is 

L m K n 20,^/1 — e\e\ 

L^7 = ^"10 (2 - e\)(l - e 2 )5/6(l - e 2 )V3 



m K n 
m H 10 

(2B\ \ 1/2 / e 4 \ ~ 1/2 e 4 

V 2 + 4(1^!) J (l-e?W_ e |)i/3-( 109 ) 

This shows that for fixed eccentricities (and fixed polytropic 
index in equ. (|109|> ). the amount of angular momentum (in 
units of Lqm) of haloes depends only on the BEC-CDM 
particle mass (in units of ma)- That means that for a given 
shape of the halo, the amount of angular momentum neces- 
sary to sustain a vortex depends on the particle mass, such 
that only above a critical mass has the system enough angu- 
lar momentum to sustain the vortex. At higher masses, even 
more angular momentum than needed can be provided. We 
show in Fig. 3, left-hand-plot, the ratio m/m,H as a func- 
tion of A for the two figures considered in this paper, for 
respective values of Lj Lqm = 1, 10, 100. The lowest curve, 
L = Lqm, establishes a lower bound on the particle mass 
for a given halo's A-value for vortex existence. We also show 
m/rriH as a function of L/ Lqm for fixed A = (0.01, 0.05, 0.1) 
(Fig.3, right-hand-plot). For a given L/Lqm, a higher par- 
ticle mass for vortex formation is required for haloes with 
smaller spin-parameter. 

However, no bounds on the other BEC particle parameter, 
the self-interaction strength g/gH, can be determined from 
these arguments, and only the energy analysis of the next 
subsection will provide us with a sufficient condition on vor- 
tex formation, constraining g/gH to be larger than a critical 
value (g/gii)crit for a given particle mass m/rriH, as we will 
show. 

We will now investigate vortex formation for the two 
models for rotating haloes described in Section 3: Halo- 
Model A is referred to a halo which has much more an- 
gular momentum than the minimum condition requires, 
L ^> Lqm- Since this is a regime which can mimic solid- 
body rotation, we will use the Maclaurin spheroid figures 
of Section 3.1. The other model, Halo-Model B, has only 
enough angular momentum to support a single vortex, i.e. 
the vortex takes up all of the angular momentum, once it 
has formed. The underlying halo will be modelled by the 
irrotational Riemann-S ellipsoids of Section 3.2 in this case. 
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Figure 3. Dimensionless BEC-CDM particle mass m/mu and halo angular momentum: Left-hand-plot: m/mu versus A-spin parameter 
for homogeneous Maclaurin spheroids (solid) and (n = l)-polytropic, irrotational Riemann-S ellipsoids (dashed) having L/Lqm = 
1,10,100 (lower to upper curves). Right-hand-plot: m/mu versus L/Lqm for A = 0.01,0.05,0.1 for Maclaurin spheroids (solid) and 
Riemann-S ellipsoids (dashed). 



4.3 Halo-Model A: L > L QM 

4- 3.1 Energy splitting and vortex ansatz for the 
Maclaurin spheroid 

For brevity, we will omit the primes on variables in the forth- 
coming sections except for the phase functions and ener- 
gies in order to avoid a confusion between rotating and rest 
frames. In order to pursue the energy analysis, we will de- 
compose the halo wave function ip into a vortex-free part, 
tpo = /e lS °, and a vortex part, w = |w|e lSl , according to 



js' 



ipow = f\w\e 



(110) 



where the corresponding amplitude \ip\ = f\w\ is a product 
state of unperturbed (vortex-free) and perturbed (vortex- 
carrying) parts, while the contributions in the phase are ad- 
ditive, S' — S'o + S[. We will insert this ansatz for the per- 
turbed dark matter halo into the energy functional in equ. 
(|107l) . This will lead to a convenient splitting of the energy 
contributions, allowing us to compare them more easily. As 
a prerequisite, we note that for vortex- free haloes, equ. (|10[l 
and (|ll|l in the rotating frame are 



A/ - f(VS' a ) 2 - f?/(m<& + gf - /*)+ 



2m 



(111) 







(112) 



+ -^-/ V5, o • (O X r) = 
and 

V- [f (vSi^fixr)] 

where is the associated chemical potential to the normal- 
ization f v f 2 = N = J v \ip\ 2 . The corresponding energy 
functional for unperturbed haloes is given by 



2m. 



(113) 



We refer the reader to Appendix B for the derivation of the 
energy splitting. It is shown there that for ip in (|110l) . the 
associated Gross-Pitaevskii energy can be written as 



£'W=£'lfe iS °}+g' f [w]-n' f [w] 
with £'[fe lS «] in lfTT3]) . 



g',[w] = I ( 2^/ 2 |V»| 2 + f / 4 (1 - M 2 ) 2 I d 3 r+ 



(114) 



+ 



and 



m 2/h . m .2, ,2 



|™| 2 [-l»-2<I>o]) 



m 



if w*vw ■ (ys' 



— ft x r d r. 



(115) 



(116) 



The terms apart from £' [fe' s °] describe the contribution of 



vortices to the energy. So, using the decomposition (|110|l 
vortices of tp (if present) are vortices of w and they are de- 
scribed via the energy functional s Q' f \w] — !Z' f \ 



similar splitting was deduced in lAftalion fc Du 



iRindler-Dalier] (|2008l ) for laboratory BECs in the Thomas- 
Fermi regime R/£ 3> 1. However, it can be easily seen that 
this splitting holds for any values of R/£ ^ 1. There is, how- 
ever, a notable change in the form of an additional term 
in <5}[ia>] which stems from the fact that the external trap 
for atomic gases is replaced here by the gravitational po- 
tential, which depends on the density profile via Poisson's 
equation and is hence affected by the presence of vortices. 
For laboratory BECs, on the other hand, the trap potential 
is fixed from the outset by the adopted laser configuration, 
and <5/[u>] does not contain the potential explicitly. 

We want to determine the critical angular velocity and 
energy above which at least one vortex in the halo will 
form. In order to derive an analytical result and to be as 
general as possible in the same time, we will consider a d- 
quantized straight, axisymmetric vortex, along the rotation- 
axis with core radius s. This core radius is given by that 
(cylindrical) radius, where the density recovers back from 
the inner vortex profile to its unperturbed bulk value. The 
most general ansatz for the wavefunction of such a vortex is 
w — \w\(r)e' d ^ , where the modulus depends only on the ra- 
dial variable. Inserting this ansatz into the Gross-Pitaevskii 
equation, one may solve numerically for the density pro- 
file of the halo in the presence of a vortex. However, it can 
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be shown that for gravitational potentials <J> falling off as 
-l/r b , b < 2 at innnitjQ this profile goes like r |d| for r -> 0, 
after imposing the constraint that it approaches the unper- 
turbed density for r — >■ oo. This is the well-known behaviour 
valid also for trap potentials used in atomic BEC gases. So, 
the density of the halo will tend to zero towards the axis of 
the vortex. For our purposes, it will be sufficient to catch 
this behaviour by the simpler profile we are going to adopt. 
Our ansatz for the vortex is 



w(r, <j>) — \w\ (r)e tSl 
with amplitude 



1*1(0 = 



for r ^ s 
otherwise 



(117) 



(118) 



and phase S{ = d<p. This amplitude reflects the fact that 
outside of the vortex, the halo (number) density is given by 
the unperturbed profile |/j 2 = po/m. Since w is going to be 
multiplied by the unperturbed wave function according to 
the above decomposition (|110[) . ip — ipow, its amplitude is 
dimensionless and the z-dependence is trivial for the straight 
vortex. We think that the consideration of a bent vortex will 
not change our conclusions given the orders of magnitudes 
involved in the final results. The constant C n is determined 



by the normalization J \ 



N. Although the mass is con- 



served, the above ansatz for the vortex causes an unnatural 
steepening or 'overshooting' of the profile at the core radius, 
which would cause a singularity in the dynamical equations. 
However, this feature is invisible to the energy calculation 
we are going to perform. 

The vortex changes the gravitational potential of the 
smooth, vortex- free halo ([62} . Owing to the linearity of the 
Poisson equation, we may decompose the total halo potential 
and mass density into unperturbed and perturbed parts, 



$ = $o + $1, 



: Po + Pi, 



where $i is the associated potential to the perturbation of 
the density due to the vortex, 



pi = p - po = Po(\w\ - 1) = 



outside the vortex 



— 1 1 < otherwise 



(119) 



according to (|117)l , and A<E> = A$ + A<l>i = 4nG(po + pi)- 
So, we are left to solve for the unknown potential of the 
vortex configuration, 

A<l>i = 4tvGpi = 

outside the vortex 

(120) 



4nGp [Cl (J) - 1) < otherwise 



The first case amounts simply to solving the Laplace equa- 
tion, while the second case is an inhomogeneous extension. 
Since we only encounter axisymmetric configurations, both 
cases can be solved analytically. However, since the solu- 
tion for d > 1 happens to be a very cumbersome expres- 
sion involving numerous hypergeometric functions, we will 
in the forthcoming analysis restrict our attention to singly- 
quantized vortices having d — 1, whose corresponding profile 



12 More precisely, we require that if 0(p — $) ~ i/r b , then < 



sC 2. 



can be found in Fig(4] In fact, physical reasoning makes the 
(d = l)-vortex more interesting, since multiply-quantized 
vortices are generally subject to splitting into several singly- 
quantized vortices. The above differential equations (|120|l 
can be solved in a standard way, and the subtle issue re- 
maining is the choice of suitable boundary conditions. For 
the outer-vortex solution, we impose that the potential ap- 
proaches a point-mass potential for large r at fixed z, or 
large z at fixed r, respectively, resulting into 

„2 



4>i o) (r,£) 



2-rrGpoCS 



(121) 



\/r 2 + z 2 

For the inner- vortex solution, we require the gradient of the 
potential to vanish, i.e. no net gravitational force, at the 
centre, resulting in 



$i\r, z) = -nGpor 2 



3©' 



(122) 



The total halo potential in the presence of this vortex is then 
given by 



4> 



(°) 

1 

(i) 



outside the vortex 
otherwise 



(123) 



with the potentials given in ([52")) . (|121l) and (|122|) . 

In Halo-Model A, the vortex is essentially considered to 
be a perturbation of the total angular momentum L of the 
halo. Since the latter is conserved, that part carried by the 
vortex, Lqm, is assumed to be small compared to L, such 
that the total angular momentum of the system before and 
after vortex formation can be given by L in (1681) . 

4-3.2 Vortex energy and critical angular velocity 

Our prescription for the halo wave function is thus tjj = ipow 
with the unperturbed wave function tpo = fe lS ° having the 
above described geometry and rotational properties of the 
Maclaurin spheroid and the vortex ansatz w from (| 11T[) . The 
normalization constant is then C„ = \J d + 1. Inserting this 
wave function ansatz into equ. (|114p will give us the critical 
angular velocity Sl c above which the energy is lowered by 
the presence of that vortex. Since our ansatz leads to an 
energy greater than that which would result if the 'real' wave 
function were used in place of the ansatz, this means vortex 
creation is energetically favoured, in general, if > Q c . 
More precisely, we have the upper bound 

JS'nl 



= £'[fe i3 '°]+ I (^/ 2 |V»| 2 + ^/ 4 (l-|» l2 ' J 



£'[fe b «]+g' f [w]-n' f [w] 



+ 



-f<fo(l-|»| 2 ) 



+ 



— f 2 I iw'Vw ■ 
m 



po\w\ 



VS' 



— n x r] £Tr. 

h 



(124) 



Since /j 2 = po/m, is constant, the integrals can be per- 
formed straightforwardly. Furthermore, VSq = since the 
spheroid has no net velocity in the rotating frame, v' = 0. 
Let us consider each term separately, where we split the 
axisymmetric integration domain according to our ansatz: 
the quantum-kinetic term is affected by the presence of the 
vortex, so 
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Figure 4. (Number) density profile of the vortex ansatz of equ, llll8l l for d = 1. 
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|Vw| rdrdzd(j>+ 



2— — c J J r = s 
2tt r s 
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(125) 



where the leading logarithmic term stems from the angu- 
lar kinetic energy of the vortex. The term due to the self- 
interaction amounts to 

' " C Jo |Ai-M a ) a «*rd0«fe = s/W ^TT (126) 



while the rotation term becomes 



K' f [w] = hf 2 



c i>1tt />awl 



c J •/ 



iw*X7w ■ (SI x r)rdrd(f>dz 



— ■Rhf 2 ca 2 dQ.. 
3 J 



(127) 



Inserting (|121[) - (|123|) and ((62j) into the second integral of 
(|115p . the gravitational potential energy due to the (d = 1)- 
vortex perturbation is finally given by 



c /i27r pa\ 1 



cJO 

c r2ir />s 



!j<f>[ o) rdrd</>dz+ 



+ 



+ 



c JO JO 



+ 



rdrd(f>dz 



= 2% 2 Gpocs 2 



2 ac . / \ 2/, / s \ 

c H arcsm(e) — c y 1 + ^-J — 



.2L^» fl *c(A 1 (e)+5/2).(128) 



The above terms constitute the total energy £'[il>] of our 
ansatz for the wave function ip = xpow. Before we collect 
the above terms, we set d = 1 everywhere and rewrite the 
semi-axes a and c in terms of the mean radius R and the 



eccentricity e. Also, we use p0[) in order to finally write the 
vortex energy in the rotating frame as 

5E 1 = g' f [w]-n' f [w] 

or 

5E' 



QqmLqm 



= 1(1- e 2 )^ 



R 



In +111(2(1-0 



13 
+ 12 



+ 



2 V n 



IQM 



(l-e 2 ) 1/3 x 



(1 _ e 2 )2 /3 + a£csm(e) (1 _ e2)1/6 _ 



(l_eY/3^ 1+ ^|y (1 _ e2)2 /3 



Hf (1 - e2)1/3 



i + ( | ) it -,-)- 



K^) 2 a) 4 ^» ,/a (-« + i 



n. 



.(129) 



For singly-quantized vortices, the core radius s is al- 
most the same size than t he healing length £ (see e.g. 
IPitaevskii fc Stringaril l|2003l )). and can be very well approx- 
imated by it, so we have replaced s by £ in the above ex- 
pressions altogether. We have also used (|39p to replace the 
coupling strength g by the healing length f. 
We can rewrite the variables appearing in (1129f) in terms of 
the BEC particle parameters m and g using (|28|l - (|3ip and 
(1406 - (14111 in order to arrive at 



5E' 



QqmLqm 



(l-O 



2x1/3 



In 



2 VW 
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+ In 2(1 - e 



a-v-i/6 
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13 
12 



3 / jm_ \ gn , 
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+ _ ^ (1 _ e2) l/3 

2 \m H J g 
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— ) (l-e 2 ) 1/3 ./l + ^(l^e2)2/3- 
9 / X 9h 
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fin 



In 



1 ( JJ}_\ ( 9h_ 
8 \m H ) \ g 



2,1/3 



1 + _y_(l_ e 2)2/3 
9H 



^i(e) + ^)-fi(e)^-(130) 



1-e 



with Q in (|65|) . Now, the total energy (| 124f) will be lower 
than the vortex-free, unperturbed energy £'[/e I,St) ] if the 
vortex energy in the rotating frame SE' becomes smaller 
than zero, i.e. if the system with vortex is energetically 
favoured. The critical condition amounts to setting SE' — 
in (|130[) and finding the respective relationship between the 



particle parameters y = rajmu, x = \J gj gu for fixed halo 
eccentricity. We see that the critical curves are just the two 
solution branches of a quadratic equation 

2/1,2(1) = b/(2a(x))[l ± ^l-Aa(x)c(x)/b 2 ] 

with 



(131) 
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,2\l/3 



-2 In x{\ 



2x2/3 arcsm(e) ,„ 
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» 2 W 6 _ 



e 2 ) 1/3 + 



l + a; 2 (l-e 2 ) 2 / : 



Ai(e) + 5/2 1 
12 

S ffi(l- 



ar* 

2x-l/3 



1/6 



12 



(132) 
(133) 

(134) 



6; 
and 

c(x) = lna; + ln^2(l-e 2 ) _ 

For a given eccentricity of the halo, these curves constrain 
the allowed space for the BEC-CDM particle parameters 
(m I ran ■, g I gii) for vortex formation, for which there is a 
minimum allowed value of g/gn > 1, according to (|40p (see 
FigEJ- This is reasonable since it makes no sense for a vor- 
tex core radius £ to be indefinitely large. Equ. (|13ip also de- 
scribes the critical curve Q — Q c in this parameter spacj^l. 

13 We can also solve for the critical angular velocity of the halo, 
fl c , above which a vortex is energetically favoured: 



n e = n QM |-(i-e 2 ) 1/3 







Hi) 


+ ln( 







+ 



+ i (1 - e2)1/3 (l) 2 (^) 2 [ (1 - e2)2/3+ 



3 + l- 



+ arCS ^ (e) (l - e 2 )V6 _ 1(1 _ e 2 )i/ 3 ^|) 2 (1 _ e 2 ) 2/ E 

-(l) 21n (f (1 - e2)1/3+ \/(f) 2(1 - e2)2/3+1 

This expression diffe rs fro m the result given in 
iRindler-Daller fc Shapirol feoich . since we have taken here 
into account the gravitational potential of the vortex more 
carefully. It can be shown that the general expression for 



(135) 



The results strongly suggest that BEC dark matter par- 
ticles with smaller coupling constant, especially those with- 
out self-interaction, g = 0, are not able to form a vortex. 

In Fig[5] (left-hand-plot), we show the relationship (| 131 1) 
(i.e. rn/m,H versus g/gn) independently of halo size, but 
for given halo eccentricity, i.e. given spin parameter, A — 
(0.01,0.05,0.1). One can easily check on which side of the 
critical curves vortex formation is favoured by determining 
the sign of SE' . The result is that for each set of solution 
branches, no vortex is allowed for parameters in the space 
outside of the region bound by the critical curves, or with 
other words, only inside that bound region, for a given A, is 
the vortex energy in the rotating frame negative. One can 
see that the parameter space of vortex existence grows as 
e or A increases, as expected. However, the energy calcula- 
tion did not incorporate virial equilibrium of the halo, which 
we have seen further constrains (m/mH,g/gH) according to 
(|81[) . We therefore also show this relationship in Fig[5]for the 
same e-values. The sensitivity to the eccentricities is weak, 
so the respective lines seem to lie on top of each other on 
the double- logarithmic plot. In light of the discussion in Sec- 
tion 2.2, we expect the virial relationship (|81[) to be valid 
only if g/gn 2> 2. Now, for a given spin-parameter, virial- 
ized BEC haloes will form vortices for those (rn/mn , g / gn)- 
values, for which the virial line lies inside the region bound 
by the critical curves. The intersection of the virial line with 
the critical energy curve defines thus a set of critical values 
(m/mH)crit, (g/gH)crit for each A, above which a vortex will 
form. The results are shown in Table [2] 
In FigO we translate the bounds on vortex formation from 
Fig[5]to the dimensional BEC particle parameters m and g. 
To this end, a halo mass and size have to be chosen, which 
will fix m,H and gn in (|48p and (|49p . respectively. We chose 
a Milky- Way halo and a dwarf-galaxy-sized halrF 4 ! For the 
latter, the parameter space of vortex formation is shifted to 
higher values of particle mass as compared to the Milky- Way 
halo. Generally, however, the parameter space of vortex for- 
mation shifts to lower values of the particle mass for haloes 
of the same size but lower mean density. 
As a further illustration of the importance of self-interaction, 
we plot the vortex energy SE' of equ. (|130p as a function of 
the spin parameter A for several chosen values of BEC mod- 
els (m/mH, g/gn) in Fig[5] right-hand-plot. Again, we see 
the trend of the result in Fig[5](lhs) at play: for a given spin- 
parameter A, the vortex is increasingly favoured for larger 
coupling strength g/gn- On the other hand, for fixed g/gn, 
a larger spin-parameter makes the vortex favoured. So, the 
total energy of the system is lowered above some critical A, 
which is smaller for higher values of g/gn- For the depicted 
examples, we see that the BEC model of the upper curve 
with (m/mH, g/ gn) ~ (100, 100) is never able to form a vor- 



cqu J135l l is a monotonically increasing function of d, so our 
restriction to the lowest critical angular velocity for which there 
is d = 1 is well-motivated, after all. In the Thomas-Fermi regime, 
as R/£ — > 00, the leading-order term in the critical angular 
velocity diverges logarithmically , f2 c ~ ln(_R/g), just as fo r 

la boratory condensates (see e.g . iLundh. Pethick fc Smith! lll997l l 
or iPitaevskii fc Stringaril 1 120031) ). 

14 Note again that the critical energy curves need R and M as 
input, while the virial constraint needs R only to specify the di- 
mensional (m, g)-space, see equ. H80H and 1821 . 
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Figure 5. BEC-CDM particle parameter space and vortex energy for Halo-Model A: Left-hand-plot: Critical curves SE' = (or 
Q = Q c ) for vortex formation in the dimcnsionlcss parameter space (g/ gjj ,m/m^) for A = 0.01 (e = 0.062) (long-dashed black), 
A = 0.05 (e = 0.302) (solid red), A = 0.1 (e = 0.550) (short-dashed blue); BEC haloes constrained by the Virial theorem: straight lines 
for the same e-values. Right-hand-plot: Vortex energy in the rotating frame in units of QqmLqm versus A-spin parameter for three 
different BEC-CDM particle models (m/m H ,g/g H ) = (100, 100), (100, 900), (100, 2500). 



tex, since it would cross the (SE' — 0) -horizontal only at a 
A-value, which is beyond the stability limit of the spheroidal 
halo, see Table I All 

We also note that the condition of a vortex to be ener- 
getically favoured, SE' < 0, does not automatically fix the 
numerical value of L/Lqm > 1- In fact, the value of L/Lqm 
varies along the critical energy curves in Fig[5] (lhs) , being 
larger for higher m/mH, according to (|108fl . The minimum 
on those curves at the critical BEC-CDM particle parame- 
ters is a factor of L/Lqm — 4 — 6 for the A- values consid- 
ered, see Table [2] This reflects the fact that the spheroidal 
haloes carry a lot of excess angular momentum over Lqm 
in all of the relevant parameter space. It seems at first 
sight counter- intuitive that (L/ Lqm)ctU decreases with in- 
creasing A. However, one must bear in mind that this re- 
sult comes from critical conditions on the energy for vortex 
formation, i.e. if A happens to be small, so must the ratio 
L/Lqm be large enough for a vortex to form, or, conversely, 
a higher A makes a vortex possible at lower L/Lqm- The 
spin-parameter, after all, does not only depend on L, but 
also on the total energy E of the system, see equ.([TJ). 

Table [2] also shows the corresponding vortex core radii 
for the chosen haloes, according to the relation in (|41[) . 
Since the vortex core shrinks with increasing self-interaction 
strength (see equ.(j39)l), the values displayed are the maxi- 
mum possible ones. 

4.4 Halo-Model B: L = L QM 

Energy and density profiles of the irrotational 
Riemann-S ellipsoid 

To account for the effects of compressibility and irrotation- 
ality, we take as our unperturbed halo with no vortex an 
(n = l)-polytropic Riemann-S ellipsoid having L — Lqm, 
just enough angular momentum for one quantum vortex. 
We will perform a similar energy analysis as in Section 4-3 
to determine if a vortex, which carries all of the angular 
momentum, is energetically favoured. Since the vortex takes 
up all of the angular momentum in this model, however, the 



Table 2. Halo-Model A: Lower bounds on BEC-CDM particle 
mass and self-interaction coupling strength for vortex formation 
in haloes of given spin-parameter; the corresponding ^-values are 
upper bounds for the vortex core radius 



Independent of halo size: 



A 


(m/m H ) crit 


(g/9H)crit 


{L/LQM)crit 


0.01 


309.41 


1.02 ■ 10 5 


5.65 


0.05 


49.52 


2549.24 


4.53 


0.10 


21.73 


454.54 


4.02 


Milky- Way-sized halo: M = 10 12 M Q , R = 100 kpc 


A 


m C rit [eV] 


dcrit [eV cm 3 


(max [kpc] 


0.01 


3.30 • 10" 23 


2.30 ■ 10" 59 


0.31 


0.05 


5.28 ■ 10" 24 


5.74 • 10" 61 


1.98 


0.10 


2.32 ■ 10~ 24 


1.02 ■ 10~ 61 


4.69 


Dwarf-; 


;alaxy-sized halo: M = 10 10 M Q , R = 10 kpc 


A 


merit [cV] 


dcrit [eV cm 3 


(max [kpc] 


0.01 


1.04 ■ 10" 21 


2.30 ■ 10~ 58 


0.03 


0.05 


1.67- lO" 22 


5.74 ■ 10" 60 


0.20 


0.10 


7.33 ■ 10~ 23 


1.02 ■ lO" 60 


0.47 



final shape will be a spherical halo with a cylindrical vortex 
in its centre. 

Before we study this in more detail, we first consider the 
energy of the vortex-free ellipsoidal halo model. The energy 
of the irrotational Riemann-S ellipsoid with (n — 1) is given 
by (|55p . where we set Kq = and use equ. (|95[l - (|98[l for the 
remaining terms in order to arrive at 

E R = Usi + W + T = 

= 4 M 2^~I^ /(ei ' e2)+ S K;(ei ' e2)Mi?2 ^' (136) 
where 
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Figure 6. Critical curves 5E' = for vortex formation in the dimensional BEC-CDM particle parameter space for Halo-Model A: particle 
mass m in eV versus g in eV cm 3 for a Milky- Way-sized halo, M = 10 12 Mq, R = 100 kpc (left-hand plot) and a dwarf-galaxy-sized halo, 
M = 10 10 M Q ,_R = 10 kpc (right-hand plot); caption as in Fig[5](lhs) 
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(137) 



is a function of the eccentricities only. Here we denned Q-g = 
y/nGp E , A = A/Qg and p B , pf" are the mean and central 
densities of the Riemann-S ellipsoid, respectively (see also 
Appendix A. 3). The total angular momentum is given in (|97[) 
and we will force it to be equal to Lqm in the calculation 
of the vortex energy in the next subsection. Although we 
will only make use of the global energies as stated in Section 
3.2 for the calculation which follows, we find it useful to 
plot the actual density profiles of the Riemann-S ellipsoids 
in the absence of a vortex, see Fig[7] left-hand-plot In fact, 
as we show in Appendix A3, analytic profiles can be derived 
if the ellipsoidal approximation is assumed. 



4-4-2 Vortex ansatz and vortex energy 

Once formed, the vortex is assumed to carry all of the angu- 
lar momentum L = Lqm of the halo. That means we antici- 
pate a transition of the Riemann-S ellipsoid to a sphere, but 
with a vortex in the centre. A spherical hal o with central 
vortex has been studied e.g. in iKain fc Lind |20ld ). where 
an approximate calculation of the density profile in the pres- 
ence of a vortex has been presented. For convenience, we will 
thus in the following take advantage of the result for the ap- 
proximate energy of a spherical halo with cylindrical vortex 
as presented ther J^l. 



15 As we indicated in Section 1, the approach in IKain fc Lind 
j201Ch to calculate the energy of the halo with vortex neglects 
angular momentum in the first place, making the very appearance 
of the vortex unmotivated. However, since we start from a halo as 
a Riemann-S ellipsoid with minimum angular momentum Lqm, 
provided by the A-spin parameter, it is meaningful to ask the 
question of whether it is energetically favoured to drive the system 
to a state having a vortex in the center. 



The usual a nsatz for a singly-qu antized, axisymmetric 
vortex is used in IKain fc Line! (20101. in our notation w = 
Iwle 1 ^. Then, the GP energy of a halo with such a vortex on 
to p of the sphe r ical b ackground density (|35p is calculated 
in IKain fc Lind (|2010t ). Using our notation, the energy of 
the spherica l halo in the presen ce of the vortex given in 
equ.(6.22) of lKain fc Linel fcoiOT l is 

Es = Usi + W + T = d-M - I ^+ 
2 2m 2 4 Ro 

+ -LQ M n Q M(l + irSi(n)(ln2 + \n(Ro/Rc)) + nT), (138) 

where Si(n) w 1.852 is the sine integral evaluated at n, 
T ~ —2.658, and pf is the central density of the halo, which 
is an (n = l)-polytropic sphere with radius Ro given in (|37p . 
Some shortcomings enter their calculation of the angular ki- 
netic part of the vortex energy: the actual vortex profile, 
as for instance given by \w\ 2 in equ. (H19[l . is replaced by a 
finitely thick, empty cylinder with radius R c , in our nota- 
tion R c — CvpVPc ■ Also, the change in the gravitational 
potential due to the vortex is neglected as can b e seen in 
the derivation of equ.(6.22) in lKain fc Linel (|2010f Fl. Using 
their result for the energy of the spherical halo in the pres- 
ence of the vortex provides us thus with a rougher upper 
bound estimate of the conditions for vortex formation, than 
if we had used a better implementation of the vortex ansatz 
of Section 4-3-1 in the Riemann-S ellipsoid. However, given 
the orders of magnitudes involved in the resulting bounds on 
the BEC particle parameters, we do not expect these effects 
to be very significant. We stress that the above expression 
for the energy, equ. (|138|) . assumes that all of the angular 
momentum is in the singly-quantized vortex, whose amount 
is Lqm = Nh. Equ. (|138p is the energy as measured in the 
rest frame, therefore E' s = Es — Q.Lqm is the energy in the 
frame rotating rigidly with Q. The energy in 136 p is also 
the one given in the rest frame, so again E' R — Er — QL is 
the energy in the rotating frame. Now, in order to compare 
the two states on an equal footing, we force the Riemann-S 



16 More prec isely, the gravitatio nal potentials Vq in equ. (6.8) 
and (6.13) of IKain fc Lind J2010T) are simply set equal, whereas 
our analysis in Section 4-3.1 clearly distinguishes between the 
unperturbed halo potential <I>o and the perturbed one, 
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ellipsoid to have the same amount of angular momentum as 
the sphere with vortex, i.e. all of the angular momentum of 
the ellipsoid shall be transferred to the vortex, leaving the 
bulk of the halo spherical. We will thus set L — Lqm in 
E' R . Moreover, since the total mass M of the halo shall not 
change across the transition to the vortex state, the mean 
and central densities, p s , p E and pf,p E , for the spherical 
and the ellipsoidal halo shape, respectively, will differ but 
be related according to 



P = P ff(ei,e 2 



.3/2 



and pf = Pcg{ei, e 2 ~ 



3/2 



(139) 



where we used p E = 3M/(4tt7? 3 ) and dToTj) (see also 
equ.(3.10) in LRS93). Furthermore, by using (fl39]i and (j36]l . 
we can write p B /pf = -^g{ei, e2) 3//2 , which we use in (I138[) 
in order to rewrite 



(140) 



\n{Ro/R c ) = \nR/i + - ln(7r 2 /3 ff ( ei , e 2 ) 1/2 ). 

Calculating the energy difference and rearranging terms, we 
get for the (dimensionless) vortex energy, 



SE' 



QqmLqm 



QqmLqm QqmLqm 



^(/(ei,e 2 )- 3/2 -l) 



J_ 

f)ii 



+ a v ; /(ei.e 2 'in -2- 
8 9h 



[<?(< 



ei, e 2 I 



'1/2 



/(ei,e 2 )] 



m 



+ ^fl(ei,e 2 ) x 



1 + ttS^tt) In 2 + TlT + E^EI m 



-^-£(ei,e 2 )p(ei,e 2 ) ( — — 
5 V m H 



3 ff (ei,e 2 ) 1 / 2 



(141) 



We stress again that formula (|141[) has been derived by 
imposing L = Lqm- In Fig(7] right-hand-plot, we show a 
spherical halo with central vortex which transitioned from a 
Riemann-S ellipsoid having A = 0.05. 

In Fig[S] left-hand-plot, we show virialized haloes which 
are irrotational Riemann-S ellipsoids according to (|103[) for 
A = (0.01, 0.05, 0.1): Thanks to the relationship in (flQ9|) , the 
constraint of having minimum angular momentum for vor- 
tex formation, L = Lqm, means to fix the value of m/mH- 
For given A and ml ran, the Virial constraint (|103j> fixes 
then the self-interaction g/gH- This means that the condi- 
tion L — Lqm can be met at only one point (for each A) in 
the BEC-CDM particle parameter space. These points are 
also depicted in Fig[S] Inserting those values for the parti- 
cle parameters in (|141|l shows that the corresponding vortex 
energy is negative, i.e. the vortex is favoured for all A-values 
of interest. Again, the vortex is favoured at high enough 
particle mass and self-interaction strength for a given A. To 
strengthen this picture by adding more 'data points', we 
plot fl 103 fl for a larger range in A according to Table lA"2l and 
the corresponding values (m/mH, g/gn) at which L = Lqm 
(see Fig(8] right-hand- plot). A vortex can be formed in all 
cases considered. In addition, FigO confirms the previous 
result of favouring vortices at high enough self-interaction 
for a given A and vice versa. Since vortex formation gen- 
erally requires L ^ Lqm, the above calculation provides 
us thus again with the critical BEC-CDM particle parame- 
ters, above which vortex formation is favoured, along with 



Table 3. Halo-Model B: Lower bounds on BEC-CDM particle 
mass and self-interaction coupling strength for vortex formation 
in haloes of given spin-parameter; the corresponding £- values are 
upper bounds for the vortex core radius 



Independent of halo size: 


A 


(m/m H )crit 


(g/gii)crit 


LI Lqm 


0.01 


44.58 


1595.07 


1 


0.05 


9.49 


68.00 


1 


0.10 


5.01 


17.20 


1 



Milky-Way-sized halo: M = 10 12 M Q , R = 100 kpc 



A 


merit [cV] 


Qcrit 


[eV cm 3 ] 


(,max [kpc] 


0.01 


4.75 • 10" 24 


3.59 


10- 6 i 


2.50 


0.05 


1.01 ■ io- 24 


1.53 


10 -62 


12.13 


0.10 


5.34 • 10~ 25 


3.87 


10 -63 


24.11 


Dwarf-galaxy-sized halo: M = 


= 10i°M Q 


, R = 10 kpc 


A 


m crit [eV] 


Qcrit 


[eV cm 3 ] 


Cmm [kpc] 


0.01 


1.50 • 10~ 22 


3.59 


IO" 60 


0.25 


0.05 


3.20 ■ 10" 23 


1.53 


10-61 


1.21 


0.10 


1.69 ■ 10~ 23 


3.87 


10 -62 


2.41 



the maximum vortex core sizes. Table [3] summarizes those 
parameters in dimensionless and physical units. Owing to 
the constraint L = Lqm, we do not calculate the critical 
condition 5E' = in Halo-Model B, since that condition is 
fulfilled at another L Lqm)- Apart from the considera- 
tion, that this model may be more realistic than the rigidly 
rotating haloes of Model A in Section 4-3, the purpose of this 
model was also to show that the results gained in Section 
4-3 are consistent with those derived here, despite the sim- 
plifying assumptions on the unperturbed (vortex- free) halo, 
on which the former were based. The same trend can be 
seen here as in Model A of Section 4-3- vortex formation 
requires a high enough mass m and positive self-interaction 
g of the BEC-CDM particles, where smaller A-spin param- 
eters require higher values for m and g. Small-mass haloes 
set tighter constraints on the particle parameters. Both halo 
models are thus qualitatively consistent with each other. 
For the same halo with fixed A, Model A, however, requires 
in general larger values of m and g then Model B. This is 
due to the fact, that Model A has more angular momentum 
L/Lqm S> 1, resulting in higher m/mH via QlOSp . which in 
turn requires higher gfgn due to the virial constraint (|8ip . 
By the same token, the maximum vortex core radii can be 
substantially larger in Halo-Model B, since the critical self- 
interaction strengths for vortex formation are smaller than 
for Halo- Model A. 

We also see from Fig[H]that a significant portion of the 
particle parameter space remains, where vortices will not 
form. In fact, for all values of (g/gH ,m/mii) which lie left 
to the depicted points on a given virial curve, depending on 
A, haloes fulfill L < Lqm- For coupling strengths which are 
high enough for them to be still in the polytropic Thomas- 
Fermi regime, g/gH 2, those haloes will just remain com- 
pressible, irrotational Riemann-S ellipsoids, and will not un- 
dergo vortex formation. 
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% =(0,0.01 ,0.05,0.1) 



1 2 3 4 5 6 




0.5 1 1.5 2 2.5 3 3.5 

r/R char 



Figure 7. Density profiles for Halo-Model B: Left-hand-plot: Profiles of the (vortex-free) (n = l)-polytropic Riemann-S ellipsoidal haloes 
having A = (0.01,0.05,0.1), employing the ellipsoidal approximation, according to cqu. llA22l l. The solid curve is the (n = l)-polytropic 
sphere, cqu. (1351 ). and is added for comparison. The densities are all normalized to = 1. The locus of the outer surface where the 
density vanishes increases with A; Right-hand-plot: Profile of the spherical halo with vortex in the centre (solid) using equ.(4.16) of Kain 
& Ling (2010) with ©o = 1 an d Rc/£, for A = 0.05. The unperturbed sphere (dashed) is added for comparison. 




Figure 8. BEC-CDM particle parameter space for Halo-Model B (no logarithmic scales here !): Left-hand-plot: Virialized haloes according 
to I I103H for A = (0.01,0.05,0.1) (lower to upper curves) with depicted points {m/ran , g / gn), according to Table[3] at which L = Lqm 
and corresponding vortex energy 5E'\ Right-hand-plot: the same but including more curves in A according to Table lA"2l and zoomed-in 
closer to the origin; dots denote BEC haloes having L = Lqm- 




Figure 9. Halo-Model B: Virialized BEC haloes according to II103D for which L = Lqm: Left-hand plot: spin parameter vs. self-interaction 
coupling strength; Right-hand-plot: Vortex energy d 141 1 in units of SIqmIjQM vs. self-interaction strength. 
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5 CONCLUSIONS AND DISCUSSION 

We have studied in this paper the conditions of vor- 
tex formation in galactic haloes composed of dark mat- 
ter particles which are able to form a Bose-Einstein con- 
densate. Haloes can then be described as fluids, obeying 
quantum-mechanical fluid equations. Quantum-mechanical 
effects make this form of dark matter behave differently from 
standard CDM, resulting in new effects with potentially ob- 
servable consequences. There are essentially two limiting 
cases one may consider. First, for quantum-coherence to 
be relevant on the scale of a halo of radius R, we must 
either require the particle de-Broglie wavelength (|21[) to 
be of the order of the halo size, XdeB ^ R, or else re- 
quire XdeB <C R but with a strong enough repulsive self- 
interaction to hold the halo up against gravity. In the former 
case, if v ~ v V i r for the halo, this translates into a condi- 
tion for the dark matter particle mass to be small enough, 
m>m H = 1.066 ■ 10~ 25 (ii/100 kpc)- 1/2 (M/10 12 M )- 1/2 
eV. In the latter case, in order for the repulsive self- 
interaction pressure force greatly to exceed that due to quan- 
tum pressure (i.e. the Thomas- Fermi regime), we must re- 
quire g > g H = 2.252 ■ 1(T 64 (7?/100 kpc)(M/10 12 Mq)" 1 
eV cm 3 . BEC haloes can then be approximated by poly- 
tropes with index n — 1. If we take R to be the radius of the 
virialized object supported against gravity by the dominant, 
repulsive self-interaction, this condition becomes a condition 
on the particle mass given by m > \/l5g/gH- 

These results apply generally to the global halo struc- 
ture. However, rotating BEC haloes add new phenomenol- 
ogy, and the possibility to distinguish this form of dark mat- 
ter from other candidates. To this aim, we have studied here 
the question of whether an angular velocity sufficient to cre- 
ate vortices occurs in BEC-CDM cosmologies. As quantum 
fluid systems, BEC haloes can be modelled as uniformly 
rotating ellipsoids, with and without internal motions su- 
perposed. To this aim, we have derived equations which re- 
late the eccentricities of haloes to their A-spin parameter. 
Once the latter is fixed, the eccentricities can be uniquely 
determined. Then, we have analytically studied necessary 
and sufficient conditions for vortex formation. Vortex for- 
mation requires as a necessary minimum condition that the 
halo angular momentum fulfills L Lqm = Nh, which 
implies a lower bound on m/m,H, i.e. on the dark matter 
particle mass. However, a sufficient condition for vortex for- 
mation can be established by an energy analysis, which aims 
to find the conditions when a vortex becomes energetically 
favoured. This results, in addition, to lower bounds on the 
positive self-interaction coupling strength of the dark mat- 
ter particle, g/gn- While the sufficient condition also re- 
quires L ^ Lqm, the amount necessary is determined by 
the model. 

We have studied two classes of models for rotating 
haloes in order to analyze stability with respect to vortex 
formation in two limits, that for L/Lqm 3> 1 {Halo-Model 
A) and for L/Lqm = 1 [Halo-Model B), respectively. In 
Halo-Model A, haloes are modelled as homogeneous Maclau- 
rin spheroids. The minimum angular momenta for vortex 
formation are then (L/LQM)crit = (5.65,4.53,4.02) for 
A = (0.01, 0.05, 0.1 ), respectively, which correspond to a con- 
straint on the particle mass m/mn ^ (m/mH)crit, where 
(m/m H )crit = (309.41,49.52,21.73), respectively. As long 



as m/mn satisfies this condition, the strength of the self- 
interaction must also then satisfy the condition that g/gn 
(<?/<m)erit, where (g/g H )crit = (1-02 ■ 10 5 , 2549.24, 454.54) 
for these same A-values, respectively. 

However, for Halo-Model B, which is an (n = 1)- 
polytropic Riemann-S ellipsoid, strictly irrotational prior to 
vortex formation, even L/Lqm = 1 can be sufficient for 
vortex formation, if the self-interaction strength is large 
enough. The condition L/Lqm = 1 fixes the value of 
m/mH for each A according to (|109l) . and the condition 
of Virial equilibrium (|103l) thereby also fixes g/gn- For 
A = (0.01,0.05,0.1), these values are given by m/mn = 
(44.58,9.49,5.01) and g/g H = (1595.07,68.00,17.20), re- 
spectively. According to equ. (|141[) . Halo-Model B makes 
vortex formation energetically favourable for those values 
of m/mH and g/gn- We interpret this to mean that, for 
L/Lqm > 1 (i.e. m/mn > {m/mn){L = Lqm)), vortex for- 
mation will also be favoured, as long as g/gn > (<?/ ' 9h){L = 
Lqm)- Furthermore, any values of m/mn and g/gn which 
satisfy the condition for vortex formation in Halo-Model 
A will automatically satisfy that found by Halo-Model B, 
which is less stringent, although more accurate. 

We can thus imagine vortex formation in BEC haloes 
composed of repulsively interacting particles as follows: If 
the angular momentum of a rotating BEC halo fulfills L < 
Lqm (i.e. if {m/mn) < (m/mii)cru for a given A, according 
to Halo-Model B), no vortex will form, and the halo can be 
modelled by a mildly compressible, irrotational Riemann- 
S ellipsoid, which has a poly tropic index of n — 1. For 
L — Lqm (i.e. if m/mn = (m/mH)crit), the irrotational 
Riemann-S ellipsoidal halo can make a transition to a non- 
rotating, spherical halo with a vortex at the center if the self- 
interaction is strong enough (i.e. g/gH = (g/gn)cru)- For a 
range of angular momenta fulfilling Lqm < L ^ 2Lqm, 
we may still expect a central vortex but now with the ex- 
cess angular momentum deforming the halo such that again 
a Riemann-S ellipsoid forms. Finally, if L ^> Lqm, oblate 
haloes described as Maclaurin spheroids can have a central 
vortex if m/m H > (m/m fl ) cr , t and g/g H > {g/gii)crit with 
the critical values now given by Halo-Model A. Those critical 
values determine as of when a single vortex is energetically 
favoured, but since L/Lqm 3> 1 [i.e. (L per particle) ^> h], 
it is also possible that multiple vortices will form, even a lat- 
tice of vorticeJH. For a dense lattice of vortices, the halo's 
vorticity will approach that of a rigidly rotating body, ac- 
cording to the quantum-mechanical correspondence princi- 
ple argum ent, applied to th e analogous problem in superfluid 
helium by iFevnmanl (|l955l ). A suitable generalization of the 
homogeneous Halo-Model A of Section 3.1 to compressible 
spheroi dal haloes, according to the ellipsoidal approxima- 
tion of lLai. Rasio fc Shapirol (|l993T ). will then constitute a 
viable model for this high-angular momentum regime. 

Generally, we have shown that BEC-CDM haloes in the 
polytropic Thomas-Fermi regime will typically form vor- 
tices, since this regime requires m 3> mu and g 3> gu, 
which largely overlaps the region of parameter space for 



It is known that laboratory BEC quantum gases confined by 
a wide range of trap potentials favour multiple vortices that are 
singly- quantized over a singl e vortex that is multip ly-quantized, 
see e.g. lAftalion fc Dul feOOlft ; iRindler-Dalled [120081 ). 
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vortex formation. By comparing the characteristic param- 
eters of Section 2.3 with the constraints derived for vortex 
formation, we see that, apparently, for these values of par- 
ticle parameters, the angular momentum of CDM haloes is 
typically of the order of the minimum value required for 
quantum vortex formation. Since vortex formation happens 
then in a large part of the particle parameter space, this 
is important to take into account when BEC-CDM mod- 
els are fitted to galactic rotation curves and density pro- 
files, especially in the very centers. While it is true that the 
vortex becomes more and more favoured for increasing self- 
interaction strengths g/gH, its size, and hence its dynamical 
importance, decreases. However, at the critical values for 
vortex formation and beyond, the size of the vortex in both 
halo models is large enough to be expected to be able to 
imprint a notable effect on the central halo dynamics, which 
can be seen from Tables [2] and [3] 

The appearance of vortices in the central parts of BEC- 
CDM haloes, whose core regions are depleted of dark matter 
mass, will change the gravitational coupling to baryons, as 
compared to a smooth or nearly-smooth dark matter distri- 
bution in the halo. As a result of depleted dark matter in the 
vortex cores, the subsequent collapse of baryonic matter can 
be delayed and so can star formation. A detailed analysis is 
necessary, however, to be able to quantify those effects. 

Our results show furthermore that haloes with parti- 
cles of high enough mass to satisfy the minimum condition 
that L > Lqm are nevertheless vortex-free, unless their 
self-interaction strength is high enough. In particular, vor- 
tices will not form in BEC dark matter which has no self- 
interaction, g — 0. Axions have often been modelled without 
self-interaction, but in doing so they will not be able to form 
vortices i n galactic haloe s , cont rary to what has been pro- 
posed in ISikivie fc Yanel \200$ ). In order to have vortices 
mimic a net rotational component, a dense lattice of vor- 
tices is needed. Both halo models studied here are in princi- 
ple able to sustain these structures. However, not only is it 
required that L 3> Lqm, but so is positive self-interaction. 
On the other hand, our analysis does not simply extend to 
cases when g < 0, such that the small attractive interaction 
of axions can be taken into account. While we do still not 
expect quantum vortices to appear in this case, other defect 
structures, like bright solitons, have been shown to being 
able to form not only in laboratory BECs with attractive 
self-interaction, but al s o in some m odels of axion haloes, see 
iMielke fc Velez Perezl l|2007l . 120091 ). 

Our results also indicate that there remain notable 
regions in (m, <?)-parameter space where vortices are not 
favoured. In these cases, rotating haloes may be modelled 
as irrotational, compressible (n = l)-polytropic Riemann-S 
ellipsoids, as we have shown. 

We have considered here the case of BEC-CDM which 
is a pure, i.e. zero-temperature, condensate. In the Thomas- 
Fermi regime, the characteristic size of virialized objects is 
in that case the radius of an (n — l)-polytrope, fixed by the 
ratio g/m 2 of particle parameters. As mentioned in Section 
2.2, this can be interpreted either as the full size of a halo 
or as the size of the core r egion of a larger halo. Recently, 
ISlepian fc Goodman! (|201ll ) have replaced the assumption 
of zero-temperature condensate by that of thermodynamic 
equilibrium in isothermal haloes at the virial temperature. 
The haloes which result contain an (n = l)-polytropic 



core of condensate surrounded by an 'atmosphere' of non- 
condensed bosons. They claim that such haloes are in- 
compatible with astronomical constraints. However, as they 
point out, their assumption of thermodynamic equilibrium 
at finite temperature breaks down if the 2-body elastic scat- 
tering collision time exceeds a Hubble time, depending upon 
the ratio a s /m with the scattering cross section in equ.©. In 
order to guarantee local thermodynamic equilibrium, they 
must assume o s /m ~ 1 cm 2 /g. In fact, this condition is 
not met by the ultra-light bosons considered here, for which 
a s /m is much smaller. Using equ.© and along with 
((57)) and (gSJ , we can write 

^ = 2.094 • lO" 102 ( — V x 
m \m H J 



R 



100 kpc 



3/2 



M 



10 12 Mi 



-5/2 



(142) 



which is valid as long as m/mn 3> 1. Even for large m/rriH, 
and for that matter, also for the critical values for the dark 
matter particle mass for vortex formation calculated above, 
the value of a B /m is much, much smaller than 1 cm 2 /g. 
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APPENDIX A: APPROXIMATE EQUILIBRIUM 
ROTATING FIGURES AS HALO MODELS 

Al General relationships 

In this appendix, we summarize some of the properties 
of the homogeneous Maclaurin spheroid and irrotational, 
(n — l)-polytropic Riemann-S ellipsoid used in this pa- 
per. A thorough treatment o f these figures of rot ation 
can be found in the works of IChandrasekharl (|l969h and 
lLai. Rasio fc Shapiro! (|l993l ) (LRS93). We make use of uni- 
formly rotating figures as models for haloes with non- 
vanishing angular momentum. We assume haloes to ro- 
tate about their z-axis with angular velocity fl = (0, 0, fi). 
Their mass density is po = m l/| 2 with / the correspond- 
ing (vortex-free) halo BEC wavefunction and their semi-axes 
(0,1,0,2,0,3) are along (x,y,z). The mean radius of the ellip- 
soids is defined as R = (aia2fl3) 1//3 , which differs from the 
equilibrium radius Ro of the spherical polytrope according 
to equ.j77j)-(l78j) (or (Kty ) for the Maclaurin spheroid, and 
equ. flTOTj) & (J37J) (or (|A13|) ) for the irrotational Riemann-S 
ellipsoid, respectively. 

The gravitational potential inside a homogeneous ellip- 
soidal body is given by 
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= nGpoiArx 2 + A 2 y 2 + A 3 z 2 - Ami - Mai - A 3 a 2 3 ) (Al) 

(see e.g. I Chandrasekharl |l969)), where the functions 
Ai,A2,Az, depending on the axis ratios 02/01, 03/01 or ec- 
centricities 

ei = sjl - (o 2 /ai) 2 , e 2 = \J\ - (a 3 /ai) 2 , 
respectively, are given by 

,oa 03 F{6,cj>)-E{e,cj>) 



A 1 = 2- 



A 2 



a 2 a3 E(8,<p)-F(e,<p)co S 2 



01 01 



A 3 = 2 



o 2 a3 gfsin0-g(g, 



(A2) 



(A3) 



(A4) 



with cos<?!> = 03/01, sin# = y fz^y^ja and the standard 
incomplete elliptic integrals 

E(6,(j>)= f (1-sin 2 9 sin 2 12 'dcf>' ', 



sin 2 sin 2 </>') 1/,2 e 



F(0,0 = / (1 
Jo 

The functions Ai,A 2 , A3 are generally defined this way, how- 
ever, equ. (|Al[ ) is only valid for homogeneous bodies. While 
we use $0 for the homogeneous Maclaurin spheroid explic- 
itly in the calculations for Halo-Model A, we do not need to 
use the gravitational potential of the compressible Riemann- 
S ellipsoid for Halo-Model B, since we make use of the 
global energy quantities, like for instance W , as provided 
byLRS93. 



A2 Homogeneous Maclaurin spheroids 

For the Maclaurin spheroid of Section 3.1, a% = 02 = a, «3 = 
c and (|A1[) reduces to (|62[) . Most of the needed relationships 
can be already found in the main text. Similarly to LRS93, 
we can write the dimensionless forms of W (|66[) and R ()77[) 
as well as L (I68[) as functions of the eccentricity only: 



\W\ 



\w\ 



3 (GM 2 /R) 



= (l-e 



2 i /e aresm(e) 



R = A = 
" R V3A 3 (e)(l-e 2 ) 2 / 3 



1/2 



and 



GM 3 R 25 



(1-e 2 ) 



-2/3^2 



(A5) 
(A6) 

(A7) 



with Q, in ([65}. On the other hand, we expressed also the A- 
spin parameter of haloes |l} as a function of the eccentricity 
only in Section 3.1.2, equ. (|84[) . Thus, fixing A determines 
the geometric and energetic quantities of the spheroidal 
halo unambiguously and independently of the BEC parti- 
cle parameters. Table I All summarizes them for a range of 
A-values which encompass those relevant for CDM haloes. 
Also, Fig lAll shows two examples of halo shapes for differ- 
ent spin-parameters. 



A3 Compressible, irrotational Riemann-S 
ellipsoids 

The Riemann-S ellipsoid rotates rigidly with angular veloc- 
ity ft — (0, 0, Q) like the Maclaurin spheroid. However, on 
this background rigid rotation, an internal velocity field is 
superposed specified by two requirements: (i) it shall have 
a uniform vorticity parallel to ft ; (ii) it shall leave the ellip- 
soidal figure unchanged, that is, the velocity vector at any 
point in the fluid shall be tangent to the isodensity surface 
passing through that point. The fluid circulation along the 
equator is given by 



V ■ dl = 7r(2 + /ij)0i02f2 



(A8) 



with fa = C'/fi. The flow is irrotational in the rest frame 
only for fn = —2. 

LRS93 have extended the analysis of I Chandrasekharl 

( 1969) from incompressible bodies to compressible ones by 
exploiting the so-called ellipsoidal approximation which as- 
sumes the following: i) the surfaces of constant density are 
self-similar ellipsoids (i.e. the axis ratios 02/ai and 03/01 
are the same for all interior isodensity surfaces), ii) the den- 
sity profile p E {m), with m denoting here the mass interior 
to an isodensity surface, is assumed to be identical to that 
of a spherical polytrope of same K p and n, but with ra- 
dius R = (ai0203) 1 ^ 3 . Both assumptions are strictly valid 
only in the incompressible limit (n = 0), but in the general 
case, where n 7^ 0, it provides an approximation to the true 
equilibrium solution. 

Specifying the spin-parameter A, we can determine the 
axis ratios (or eccentricities, respectively) using 106 p and 
equation (5.16) in LRS93, which is 



4(o 2 /oi) 2 



(1 + (02/01)2)2 

4Bi 2 (o 2 /oi) 2 



(o 3 /oi)2A 3 - (02/01)2 ^;-^^ 



1 + (o 2 /oi) 



+1 = 0.(A9) 



The constants K n and q n in the energy terms appearing in 
Section 3.2, depend on the polytropic index n via 



5 J X1 e n x 4 d x 
3 Y?\e>\ 



and 



= 1(1 



= 1 
= 1(1 



for n — 
0.653 for n = 1 



(A10) 



for n — 



r) w 0.523 for n = 1. 



(AH) 



6 = p s I pi and X = r/R char , with R char = s /K p /{2ttG), 
are the dimensionless Lane-Emden variables for the density 
and radius, respectively, of a polytrope. We denote the first 
zero of the density profile as xi 0- e - Xi = Ro/ Rchar) and 
0i=O(xi)- 

Now, specializing to n = 1, the dimensionless versions 
of W, R, L in (|95[) . (|101l) . (|97[) can again be given as functions 
of the eccentricities only: 



I IT I 



I Ml 



3 (GM 2 /R) 
with /(ei,e 2 ) in jSBJ 
R 



= /(ei,e 2 ) 



R 



Ro 



[/(ei,ea)(l-2t)] 



-1/2 



(A12) 



(A13) 
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Table Al. Parameters of the homogeneous Maclaurin spheroid as a function of A 



A 






t 


fl 




L 2 




IV 


H 





01 


,6240952030(-l) 


.5204(-3) 


.4558481107(- 


1) 


.2500065974(- 


3) 


.9999996616 


1 


000520746 





02 


.1242809380 


.207472(-2) 


.9086159255(- 


1) 


-1001034014(- 


2) 


.9999946162 


1 


002083877 





03 


.1850994726 


.464440(-2) 


.1355078928 




.2255299174(- 


2) 


.9999730069 


1 


004690543 





04 


.2443945077 


.819899(-2) 


.1792395408 




.4016798968(- 


2) 


.9999158013 


1 


008343664 





05 


.3017559569 


.1269872(-1) 


.2217968643 




.6291209890(- 


2) 


.9997977879 


1 


013048263 





06 


.3568449744 


.1809622(-1) 


.2629514670 




.9085911668(- 


2) 


.9995887660 


1 


018812234 





07 


.4093987579 


.24338998(-l) 


.3025107370 




.1241015884(- 


1) 


.9992548016 


1 


025647448 





08 


.4592299927 


.31371503(-1) 


.3403173945 




.1627516169(- 


1) 


.9987594192 


1 


033570749 





09 


.5062220297 


.39137176(-1) 


.3762479967 




.2069422508(~ 


1) 


.9980646348 


1 


042605066 





10 


.5503211259 


.47580078(-l) 


.4102099657 




.2568290981(- 


1) 


.9971317748 


1 


052780577 





20 


.8470361891 


.1590675677 


.6351541756 




.1124337842 




.9662629271 


1 


231978660 





30 


.9658444276 


.3004590826 


.6495893878 




.3065151508 




.8638313134 


1 


703156470 




(see HDII), and 



L 



3 ( K\ 



(y) 7 2 (ei,e 2 )fi 2 



GM 3 R 4 
where 

7(ei,e 2 ) = 

\n 2n-1/3/, 2-.-1/3 . /-, 2x2/3/! 2s- 

[(1-ei) (l-e 2 ) +(1-61)' (l-e 2 ) 



(A14) 



(l- e ?)V3 (1 _ e | } -2/3 



'(1 - e 2)-V3(i - e|)-V3 + (1 _ e f)2/3(i _ el)- 



1/3 



(A15) 



and Q, in (|94[) . The angular velocity of the internal motions 
A is related to fi via (|93[) . defining A = A/fla- The A-spin 
parameter determines the above quantities unambiguously 
via equ. (|106[) . Table [A"2l summarizes them for a range of A- 
values, encompassing those relevant for CDM haloes. Note 
that A is a monotonically decreasing function of A, while Cl 
is not monotonic, in contrast to the case of the Maclaurin 
spheroid. Two illustrative examples of halo shapes can be 
found in FigfA^I 



A3. 1 Density profile for n — 1 

LRS93 confirm their results by numerical comparison of the 
density profiles of rotating versus non-rotating configura- 
tions for different polytropic indices n. However, we observe 
that the ellipsoidal approximation in conjunction with the 



known analytic density profile of the equilibrium (n = 1)- 
sphere, equ, (135[) , makes possible the derivation of an ana- 
lytic profile for a rotating, ellipsoidal (n — l)-polytrope, as 
follows. 

According to the ellipsoidal approximation, p s (m)/p s 
and p E (m)/p E have the same dependence on m/M. From 
this follows that the shapes of the cumulative mass profiles 
are the same. The ellipsoid has a larger volume, however, 
than the equilibrium sphere according to (|101[1 . Since we 
consider the case where the total mass M shall be the same, 
this means that the mean densities (and hence central den- 
sities) are different according to 1)1390 . 

Now, the equation of an ellipsoidal isodensity surface 
can be written as 



2 i 
X + 



V 



+ 



(A16) 



whose dimensionless form, after dividing by R c har, will be 
denoted carrying a 'tilde'-sign as 



q = q I 'Rchar- 

The outer surface of the ellipsoid is then given by 



(A17) 



-2 -2 

-2 _ ~2 , V , z _ "1 / « , o\ 

Umax — X T 2 + 1 2 ~ ~n2 ' l A18 J 



1-ei ' l-el Rl har 



Writing the semi-major axis in terms of mean radius and 
eccentricities equ. (|87[l . and using (|10ip . we see that 



qmax 



Tvg(e 1 ,e 2 



-1/2 



(l-eDV^i-e^i/e- 



(A19) 
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Since the isodensity surfaces are self-similar ellipsoids, it fol- 
lows more generally that 



1 = 



Xff(ei,e 2 ) 1/2 
(l-e?)V6(i_ e 2)i/6 ! 



(A20) 



with x from above. Because of ()101[) . the volumes differ ac- 
cording to V E — V s g(ei, e2)~ 3 ^ 2 and correspondingly we 
have for the volume elements that 

dV E = dV s g(e ll e 2 y 3 / 2 

= 4nRl har (l - e?) 1/2 (l - e 2 2 ) 1/2 fdq. (A21) 

Using now the fact that the shapes of the cumulative mass 
profiles are the same, we can derive the following density 
profile of the (n = l)-polytropic ellipsoids, 



E / ~\ E 
P (?) = Pc - 



a [ 9 -(l- e ?) 1 /6(i_ e |)V6 ff(ei)e2) V2] 
g(l-ef)Ve(i_ c 2)i/6 fl(eiie2) i/2 



(A22) 



with q from above and with pf , as usual, denoting the cen- 
tral density of the ellipsoid. Using this density profile, the 
total mass can be calculated to be 



M 



: 47r 2 i? c 3 w/9 f 5 ( ei ,e 2 )- 3/2 = 4it 2 R 3 char pf = M s , (A23) 



i.e. the mass of the (n = l)-polytropic ellipsoid is the same 
as that for the equilibrium (n = l)-sphere, as we have de- 
manded. 



APPENDIX B: SPLITTING OF THE 
GROSS-PITAEVSKII ENERGY FUNCTIONAL 

In this appendix, equ. (|114p to (|116[1 are derived. All energies 
and wave functions, tp, f, w, as well as the respective ampli- 
tudes and phases, S, So, S w , are those in the rotating frame 
of reference, so primes are omitted in this section altogether. 
The integrations are over the whole space. 

We insert (|110p into the energy functional (|107l) consid- 
ering each term separately: the kinetic energy term has the 
following form, 

/ £ |vv>i2 = / £{/ 2 lvH 2 + M 2 [/ 2 (vso) 2 + (v/) 2 ] 

+^V(/ 2 ) ■ V(M 2 ) + / 2 V5 ■ {iwVw* - iw'Vw)^ . (Bl) 
The rotation term is given by 
-ih J ■ (VV> x r) = 

= ih ( fw*Vw (nxr) - h [ f\w\ 2 VS ■ (fi x r), (B2) 



while the potential energies due to gravitation and self- 
interaction are simply 

/ + §!<)=/ (y<J>/ 2 H 2 + f/ 4 H 4 ). (B3) 

It is advantageous to write the energy such that the vortex- 
free contribution is clearly separated from the other terms. 
Using the above expressions, we may thus recast the func- 
tional into the following form 

fc2 



£[iP]=£[fe lS <>]+ J ^(M 2 -l) 



(V/) 2 + / 2 (VS ) 2 - ^/ 2 VS - (nxr) 



+ 



+ I ~r [*H " *o] + / 4^V(/ 2 ) ■ V(\w\ 



£/ 2 |vH 2 + f // 4 (IH 4 - ih- 



+ / tt~ [/ 2 V5 • (iwVw* - iw*Vw)+ 
2m 



i—z- f 2 w*Vw ■ (ft x r) 



(B4) 



with £ [fe tS °] being the vortex-free energy of equ. (|113|l . Eval- 
uating the second term in line 3 results in 



v(/ 2 )-v(M 2 ) = j v(/ 2 )-v(M 2 -i) = 

- j (\w\ 2 - 1)A(/ 2 ) = -2j(\w\ 2 - 1)(/A/ + (V/) 2 ), (B5) 
which, by using becomes 

/£v ( f).v(H 2 ) = /£(H 2 -i)x 

^/ 2 VSo-(fixr)-/ 2 (V5 ) 2 - 



9m 

-_/ 2 (m$ + 3 / 2 -^)-(V/) 2 



(B6) 



assuming that the halo mass density either goes to zero as 
r — ^ oo or is identically zero beyond some finite radius. 
Furthermore, we take advantage of rewriting the following 
terms: 



9 
2 

and 



/ 4 (M 4 -i) = <?|/ 4 (M 2 -i) + § |/ 4 (i-H 2 ) 2 (B7) 



2m 
h 2 



f X7So ■ (iwVw* — iw*X7w) + ihf w*S7w • ($1 x r) 



J ifw'Vw ■ (vSo - x r 



(B8) 



respectively, using the conservation of particle number, 
J (\w\ 2 — l)f 2 p = 0. Collecting all of the above expressions, 
we shall finally arrive at the following splitting of the energy 
functional: 



£lTP]=£{fe iS °} + g f [w]-n f [w], 
where 

g,[w] = [ (|-f|v w | 2 + f/ 4 (i-H 2 ) 2 ) - 



(B9) 



/ (f/ 2 $o + f/ 2 H 2 [*-2$o]) 
d 

UAw] = — [ if 2 w*Vw ■ (VSo - ?H x r 
m J ' \ a 



+ 

and 



(BIO) 



(BH) 
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Table A2. Parameters of the irrotational, (n = l)-polytropic Riemann-S ellipsoid as a function of A 



a 2 /a 1 a 3/ a l e l e 2 






01 


.7073606352 


.8191876671 


.7068528360 




5735255584 


.8990308186(-3) 





02 


.6157806530 


.7444171156 


.7879176273 




6677148778 


.3256446859(-2) 





03 


.5551042656 


.6884452492 


.8317807730 




7252883143 


.6676951963(-2) 





04 


.5095561821 


.6425906919 


.8604373872 




7662096337 


.1087731987(-1) 





05 


.4731628114 


.6033952054 


.8809750019 




7974423027 


.1565115011(-1) 





06 


.4429370293 


.5690393124 


.8965527247 




8223103191 


.2084571079(-1) 





07 


.4171522426 


.5384160866 


.9088366226 




8426791309 


.2634641496(-1) 





08 


.3947154535 


.5107881879 


.9188034125 




8597065936 


.3206621692(-1) 





09 


-3748897504 


.4856329988 


.9270694014 




8741627940 


.3793824738(-l) 





in 


.3571546372 


.4625629986 


.9340452693 




8865864156 


.4391062656(-1) 





20 


.2417900586 


.3025300680 


.9703285874 




9531398418 


.1032663262 





30 


.1766421365 


.2101855698 


.9842751422 




9776615090 


.1550787050 


A 




a 


A 


L 2 




\W\ 


R 





01 


.7169963591 


.6760714962 


.8424269299(- 


6) 


.9919850657 1.004935593 


0. 


02 


.7176346622 


.6408208912 


.3253671535(- 


5) 


.9842263532 1.011279929 


0. 


03 


.7165548594 


.6081343468 


.7088561436(- 


5) 


.9766699520 1.018697841 





04 


.7140826379 


.5777254698 


.1223320264(- 


4) 


.9692743288 1.026958151 


0. 


05 


.7104728779 


.5493488048 


,1859786170(- 


4) 


.9620070989 1.035897363 





06 


.7059271906 


.5227939592 


.2611162661(- 


4) 


.9548427521 1.045397324 


0. 


07 


.7006067868 


.4978802570 
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Figure A2. Irrotational Riemann-S ellipsoidal halo rotating about the z-axis having 123 = 1 and A = 0.05 {left-hand-plot) and A = 0.2 
{right-hand-plot), respectively. 



APPENDIX C: THE CHEMICAL POTENTIAL 
OF BEC HALOES 

We have dealt in this paper mostly with the energy of a given 
BEC halo. The equation of motion for stationary systems in 
equ. (|51[) also involve the chemical potential p, which is fixed 
by the conservation of particles. For a static ground state 
without vortex, (|51[) easily leads to the following expression 
for the equilibrium condensate chemical potential 



M = 



h 2 A^/n 
2m -v/n 



+ m$ + gn 



(CI) 



with corresponding time-independent particle number den- 



sity n = 



Multiplying (JSl]) with dtp* /dN results in 



the well-known thermodynamic relationship, p — dE/dN, 
where E is the total energy of the BEC halo wave function 
under consideration. In the TF regime, for instance, (|C1|) 
reduces to p = gp/m + m$. Multiplying by p and integrat- 
ing results in fi = 2mE/M, where E is the energy in the 
Thomas-Fermi regime, i.e. (|55[) with Kq — (and T = if 
v = 0). The chemical potential for (n = l)-polytropic BEC 
haloes is thus given by 



H = 2 



TV" 



(C2) 



We can derive this relationship also by using the energy 
expressions of LRS93. The total energy of the polytropic 
sphere with index n according to LRS93 is 

E = U + W = kiK( P c f /n M - k 2 (p! ) 1/3 GM 5/3 , (C3) 

where the constants k\,k2 depend on n, xi an d 8[. Using 
the formula for the central density of the sphere (see also 



Xi M 

47T|0i| R 3 



M 



we rewrite 

„ , v i /w {Nmf' n+1 1/3G,.. a 

The chemical potential of this sphere is thus 

dE (1 \U JV 1 ( „ U „, 
"=m={n +1 )n +2 N = N ( E+ n +W 



(C4) 



(C5) 



(C6) 



which reduces to p = 2E/N for n = 1, in accordance with 
(IC2l). 
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